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■ABSTRACT 

THE  CQOTROL  OF  HEAVE  AHP  PITCH  CF 
A  SEMISUBMARINE  IH  AS ’Em  REGULAR  WAVES 

The  purpose  of  this  investigation  was  twofold.  First,  mt'hem&tical 

expressions  were  sought  to  represent  the  coupled  heaving  and  pitching 

motions  of  a  semi  submarine  in  astern  regular  waves.  For  the  sake  of 

simplicity,  the  linear  theory  was  used;  but  corrections  were  made  for 

scsne  nonlinear  effects.  Since  a  complete  analytical  setup  of  the 

motions  did  not  give  results  that  corresponded  well  with  experimental 

values,  so'ie  correction  factors,  consistent  with  the  theory,  were 

introduced.  In  this  way  reasonable  correlation  was  achieved  in  the 

range  of  experimental  results.  Once  the  analytical  expression  of  the 

motions  was  found,  a  mathematical  model  was  set  up  with  a  control 

system;  and  the  motions  in  waves  were  simulated  on  the  computer. 

Different  types  of  controls  were  studied,  which  led  to  a  selection  of 

an  optimum  control  device. 


-  it  * 

TABU  or  COKBPBg 


ABSTRACT  . . . . . .  i 

tabls  cp  coHsnrs  . . . . .  u 

U9T  OP  FXGESS3  . . . .  iU 


TABU  OP  HCMESCIAIURS  lv 

AtSSttOEBGtSST  . . . - . .  is 

X.  XSTROSXSCTSQB  l 

XX.  BQOmCBS  OP  IK0XOR  . . k 

xxx.  ammm  «  tks  wanness  gp  kotich  ..............  23 

XV.  TfiS  CCSiTRQS.  CP  3EAVE  ASS  PlTCfi  . . SB 

v.  concLosiosa  abd  sfficttemmcm  . . 42 

AKPBKDIX  I.  OFFSETS  OP  SBW3UBMARXBS  ...............  45 

APIHBXX  XX.  TBB  IBVXBKX  CP  3USCK  00  m  B^WIXCBS 

€P  MOTXflK  . . 1,6 


AFPSHDXS  XIX.  Bccrnno  FORCES  AS®  SttKBTS  * . «...  43 

AHKHDXX  XV.  THE  Bgn<SVBU8G  EE8IVATIVBS  AW  TEE 

SOLJmOES  TO  THE  EQUATIONS  OF  M2PIOK  ....  51 

APFK8DXX  V.  EXPECT  OP  ERXCTXCB  00  THE  S8QTXCSS  OP  A 

SBMX&JBX&XHE  . 57 

ABVSBCDXX  VX.  T2S  COKCRC*.  Ot®FACSS  ...................  59 

Apwacat  vxx.  xebal  ccwerol  . . •••• .  63 

AWEHDXX  VXXX.  COBZRCL  3I8TKM3  . . . .  66 

AVRBDXX  XX.  COMPUITO  PROGRAMS  . . .  ?6 

FIGURES  . 86 

BXXJOORAHT  . 95 


ill  . 


Liar  cp  mum 


Fisoxn  Bum 

mnmmrnmmm  >1  miWliil 

1  SUmneicm  car  t»ffc  3a&l  . . . .  66 

2  Sending  CoefTidenti  fcs?  Ssada g  sad  Htcfe&ng 

in  u*we»  of  aofl  6-£t  Length  . . .  85 

3(A)  SMw  Aeplitatfifta  la  4*  sad  B~tt  timm  •«*«**•  90 

3(1)  Pitch  AoplltutiM  in  4®  asst  S-ft  tooa  .  91 

4  Urinating  Aog&ee  sad  Haemtte  ...*•.•.»»»»««.*.  22 

5(A)  aeewe  Aoqditts&ts  viifc  Qxxtoeoh$  S>£ 

5(1)  Pitch  toplltuflea  dth  CosalaroS.®  •»»••«•  .*.«•••  $4 


TABLE  Cg  HCHMCUTCRE 


Ally  eystools  used  more  than  once  in  the  text  are  included  in  this  table; 
other  syabol*  are  defined  stars  used.  In  practically  all  cases  diasmelonal 
Visntitlee  ere  lapOied  in  foot,  pound,  and  second  units. 

a  *  virtual  rasa  la  heave  aiuatico  or  wra  a^pUtua*  at  the  free  sorfai 

a  virtual  mass  or  control  fins 

a^  m  wvt  aeplltudc  at  depth  h  halcw  surface 

a^  2,  3  *  dlaumsieaa  sf  body  of  isolation 

A  «  virtual  racment  of  inertia  ia  pitch  equation 

A  (in  chapter*  on  eaustrols)  »  area  of  control  fins 

-V  *  virtual  aaaeat  of  Inertia  to  the  gressnoa  of  tha  control  folia 

A(s}  *  tmansvaree  cress- ss oti oosl  asaa  mt  distance  x  frraa  nidflhlpa 

k"  *  sercsss,  of  vaberplsa® 

b  *  heave  djsspia^  confidant 

«  heave  dating  coefficient  des  to  pseaeac®  of  foils  only 
B  «  pitch  dashing  coafflaient 

B£1  *  pitch  dsaaplug  coefficient  du®  to  prassnee  of  folia  only 
S *  »  breadth  of  strut 

e  m  buoyancy  coefficient  ia  heave  equation 
c  <*  wave  celerity 

C  »  coefficient  of  longitudinal  righting  accent 
Cjj  »  foil  drag  coefficient 

*  coefficient  of  longitudinal  righting  iscraent  due  only  to  the  proaeaee 
of  control  fins 

*  foil  lift  coefficient 

Cg  a  noollnsar  spring  coefficient  ia  pitch  equation 

&  »  coupled  virtual  saooent  cf  inertia  coefficient  in  heave  expiation 

4=  *»  coupled  virtual  moment  of  inertia  coefficient  due  only  to  the  presence 
of  the  control  foils 


-  V  - 


D  -  coupled  virtual  mass  coefficient  in  pitch  equation 

Df  <a  coupled  virtual  mass  coefficient  due  only  to  the  presence  of 
control  fin3 

e  =  coupled  heave  damping  coefficient 

ef  a  coupled  heave  damping  coefficient  due  only  to  the  presence  of 
control  fins 

E  a  coupled  pitch  damping  coefficient 

Rp  s  coupled  pitch  damping  coefficient  due  only  to  the  presence  of 
control  fins 

f,,  fg  =  coefficients  of  cos  out  and  sin  mt  in  expression  of  exciting 
force 

fH*  function  of  time  lag  in  expression  of  varying  anti-heaving  foil 
angle 

f  a  function  of  time  lag  in  expression  of  varying  antl-pitching  foil 
P  angle 

Fr  =  Froude  number 

g  ss  coupled  spring  coefficient  in  heave  equation 
g  e  gravity  acceleration 
G  s  coupled  spring  term  in  pitch  equation 
h  j=  depth  coordinate  below  surface 
I  s  longitudinal  moment  of  waterplane 

L*  Ja  maximum  moment  of  Inertia  of  semi  submarine  ax-ound  the  y-axle 

k-,  o  o  t  »  control  parameters  in  expressions  of  varying  angles  of 
*  *  *  control  fins 

kjj  •=  ratio  of  added  mass,  a,  to  actual  mass,  ra 

kp  ss  ratio  of  A  to  J 

L  =  ship  length 

m  a  mass  of  the  ship 

*1>  ®2  *  coe^iclento  of  cos  ait  and  sin  <ut  in  expression  of  exciting 
pitching  moment 

Mq  =  maximum  exciting  pitching  moment  on  semi. submarine  due  to  strut  sad 
body 


-  Vi 


Up  a  pitching  moment  produced  on  acmsioubrnsrina  through  the  control 
r  fine 
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attach 

^xcit^^pitchiisg  acxKtnt  on  submerged  body  of  revolution  due 
to  waves 


M  «  exciting  pitching  aoseat  trots  variation  ia  buoyancy  due  to  the 
*  waves 

Mo  »  exciting  pitching  moment  directly  due  to  stmt 
V  %>  \>  K>  \>  %  *  hydrodynamic  derivatives  defined  in  Ref.  21. 

JS(x)  »  heave  tearpiag  force  a?  a  strip  of  unit  length  at  distance  x 
from  midehip 

0  s  origin  of  axis  eyatea,  canter  of  body  of  revolution  (also  approximate 

position  of  center 

p  s  hydrodynamic  presets®  °*  S^svlty ) 

K  ts  r(x)  <a  radius  at  body  of  revolution  at  distance  »  tram  midship 
t  a  time  coordinate 


tj_*  tjj  »  time  lag  of  control  systsao 

=  natural  heaving  pariod 

Tp  a  natural  pitching  period 
2tt 

T  a  mm.  *  period  of  encounter 

IV 

tj  a  absolute  value  of  speed  of  wave  particle  with  respect  to  the 
control  foils 
V  a  ship  velocity 
v  «  wave  particle  velocity 


V* 


wave  particle  velocity  at  depth  is  + 


l. 

10 


v  0«  wavs  particle  velocity  at  depth  b 
w- 

W  a  Z  *  h8SV8  velocity 


w(x  )  a  B*  j»  local  width  of  strut 

W  »  work  done  by  e  control  foil  during  s  quarter  cycle 

x  *  coordinate  along  longitudinal  exis  of  aemiBubniarine  with  origin 
midship 
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£  -  phase  8 ngle  between  pitch  m&  the  wr.va  motion 

m  phase  angle  boivean  ©acitisg  pitching  msraent  end  veva  motion 

«B  %  piuuja  engl*  betwwta  eseitisas  fasave  force  wad  w*v®  motion 

6  a  pitch  reapane®  of  the  vessel 
0O  *  maximum  pitch  amplitude 
1  c  wssve  length 

fx  ~  phase  leg  bsteeen  the  motion  of  atiti-tof&vln®  foil  ead  the  wres 

p  »  water  density 

v  s  phase  angle  between  motion  of  aati-pltcMng  fin  em&  ttsa  waves 

p^s  »'«  potential 

g?r 

in  e  -»—»(¥  *  e)  ®  frequency  of  encounter,  cycles/sec 
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natural  frequency  of  the  aenil submarine  in  heave  or  pitch. 
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a.  nfg-.ooucTioa 


She  sc- Gelled  gsadsubraarisas  (Mg.  1)  consists  of  a  fully 

eutorus5?gcd  body  of  revolution  with  a  rolabl veXy  large  bestow  strut 

whose  top  ri&as  out  of  the  water. 

Inters  at  in  this  type  of  vessel  vag  aroused  beeawe®  it  aeaawd 
some  of 

to  hol&g,the  QSveatagss  of  botfe  a  submarine  and  a  surface  ship, 

A  subtsorin©  is  sot  subject  to  ths  Lwgs  forces  m&  sjosGeats  that 

are  caused  by  &,  Cheng®  ia  buoyancy  due  .to  wave®  and  by  t&®  wavs 

orbital  velocity.  Also,  the  resieteac©  sfes  to  maressfeiag  by  tfes 

vassal  is  .less  for  a  subaisreod  body,  whils  its  increase  in 

at  high  speeds 

frictional  resistance  is  not  as  gj\a&t|*& a  tSse  fiaersaaa  in  wave 
realsfcen.ee  (Kef.  l).  On  the  other  toed,  the  surface  ship  tea 
gerraffiaeat  access  to  the  air  wfeieih  is  ®a  advantage  from  a  psycho¬ 
logic®!  viewpoint  ea&  with  raapeefe  to  issfcalletias  of  relatively 
low- coat  als>bre&tSilng  macMnsry. 

EKperiiaeafca  have  been  csrried  out  at  the  H.X.T.  and  other 

of  -  the  sesalsubmarine 

towiag  tasks  to  investigate  the  p®sr£os®aaee,.in  calm  water  and  in 
regular  waves.  3h®  motions  in  the  vertical  plane  (pitching, 
heaving,  end  surging)  received  special  attention,  l&e  teats 
(Kef.  3)  showed  that  the  serai  oubsaarged  body  behaved  very  well  in 
ahead  but  poorly  in  astern,  waves.  'Ibis  can  easily  be  understood 
since  the  natural  XrsKjuency  of  both  heave  end  pitch  are  small 
(see  App.  2)  eo  that  synchronism  with,  waves  occurs  only 
ia  astern  seas,  \ hile  in  a&aad  tseas  the  dlispsrity  between  the 


2 


natural  and  the  encounter  frequencies  becomes  larger  and  larger  with 
increasing  forward  speeds. 

It  was  tfca  author's  intention  to  find  a  theoretical  foundation 
In  expressing  these  motions  ra&thematically,  such  that  It  would 
become  possible  to  predict  the  motions  of  the  semisubmarine  in 
regular  waves  end,  by  superposition,  in  confused  sous. 

Moreover,  with  the  analytical  sathods,  it  Is  possible  to  find 
the  properties  tsaoasaaxy  to  control  excessive  motion.  Control  of 
heave  and  pitch,  with  this  type  of  vessel,  doea  not  eoera  to  be  m 
difficult  e  task  ss  it  is  with  surface  ships.  Xadaod,  as  was  mentioned 
before,  the  forces  and  actaanfcs  noting  on  the  sejaieuaeariae  will  not  be 
an  large  as  on  the  surface  ship  while  the  frequency  of  motion  at  which 
control  is  necessary  is  much  lower  than  with  surface  vessel®.  3Mo  is 
due  to  the  low  natural  frequencies  of  th®  aasaieubeasrged  ship. 

Thee?  considerations  promoted  the  desirability  of  investigating 
the  motions  analytically.  However ,  literature  on  the  motione  of  sear 
surface  vessels  ia  quite  rare  end  that  which  ie  available  (Ref.  2?) 
proved  inadequate  for  the  purposes  of  this  paper. ..  Hevsrthelasfi,  theories 
have  been  developed  in  the  past  years  for  surface  ship  and  daeply  sub¬ 
merged  submarine  motiomi  that  bear  much  relsvaace  to  the  semitebaasrinc . 

In  this  work  an  attempt  has  been  made  to  combine,  modify,  sod  adapt 
these  theories  in  such  a  way  that  an  analytical  expression  for  heaving 
and  pitching  motion  of  e  eeraiaubmarlns  is  obtained. 

He  .large  surge  motion  observed  experimentally  (Ref.  3)  must  have  a 
considerable  influence  ca  these  ship  motions  since  the  surfaos-pierclng 


-  3  - 

sail  establishes  a  relatively  long  moment  arm.  even  the  small 

surge  force  induced  by  the  vaves  on  the  strut,  or  by  the  motion  of 
the  body  itself,  must  have  considerable  effect  on  pitching  and,  by 
coupling,  on  heaving.  Hie  surge  exciting  force  due  to  the  waves  on 
the  strut  has  been  taken  into  account  in  this  work  but  not  the  effect 
of  the  self-induced  force  by  surging  on  either  the  strut  or  body.  This 
effect  would  cause  the  motion  to  become  an  involved,  nonlinear  vibration 
problem  (see  App.  II )  in  which  many  properties  have,  as  yet,  not  been 
sufficiently  Investigated. 


13»  motion*  in  the  vertical  plme  of  a  body  la  a  fluid  can  be 
exjareugod  as  a  sat  of  three  differential  equations  in  which  the 
heave,  pitch  and  marge  are  coupled  'with  each  othss?.  'M  the  earga 


ie  emitted,  two  e^patio&s  remain.  They  as?©  of  the  form: 

(*„.4>£  +  «i  +  c*+40+t&  +  40  -  “I-) 

MjvaJi?  +Sfr  +  ce+i)i,Ei/e*sfi  ■«• 


there  M 


x& 

,£■  Ze  , 


represent  the  glass©  Mgs 


with  respect  to  tfca  origin  cf  ties. 


Tb^  asie  sys to,  tMeh  moves  with 

tbs  fcsdy,  is  ri^it-tesSsd.  ra©  crigan 

ia  at  cesasteJt  dsplti  sepal  to  h;  x  i©  in 

the  direction  cf  the  longitudinal  esd,® 

of  the  body  as  revolution  at  sreet, 

ooeitiv®  to  tSsa  bow,  while  y  is  to  port  end  s  u$&&d8*,,  $& s  directions  of 
{note  that  these  directions  are  opposite  to  those  specified  in  Ref,  21). 
this  awing  coordinate  eystesn  do  not  cfosag®  duslsg  th©  motion, 

The  aquations  (l)  Si’®  2TO5i0O3i iSfe30  tbs  motions  eg®  fespt  small, 

m  waa  shown  by  $b&wwits  (Raf .  4)  who  dsvalcpsd  the  sans  elation®  from 

Hewton'ti  basis  force  equations  by  means  of  Taylor’s  expansions  {the 


symbols  are  those  of  Ref,  4) :  / 

■  (fn-ZcJz  - i(ir*  -Z**  +  VzJjB sU^6%(y) t-£z) 

\l1  -ty  j?  -(M|  +Wfji  ~(fy  +  v  -%2  -  Al^ri  -  ht  z  *  ft  cos  C tot  -  e*  ) 

wherein  w  *  *,  q  *  e,  end  ®y  «  J.  Q»  coefficients  lil»  ^  and  ^  are 
"hydrodynamic  derivative* "  (see  Ref,  4). 


*  5  “ 

Tbmat  derivatives  have  to  ba  found  togattor  with  the  magnitude 
of  the  ferae*  and  moomta  that  regular  waves  taduea  oa  tba  semi submarine 

The  Bscitlag  Forces  and  Mxacnts 

fcylov,  in  hie  treatise  on  ship  notions  around  the  turn  of  the 

century,  made  the  nov  famous  assumption  (the  so-ccUsd  Frouda-Krylov 

hypothesis)  that  the  configuration  of  the  flow  In  a  wove  le  not 

influenced  by  the  pressure  of  a  body  in  it  so  that  the  expression  for 

the  farces  end  moments  that  the  mewing  vats?  exert*  on  the  moving  chip 

is  derived  primarily  from  buoyancy  consideration*.  However,  this 

assumption  i#  not  entirely  right  since  it  is  known  frea  experiments 

mod  theoretical  hydrodynamics  that  a  body  moving  through  a  fluid 

change*  tbs  potential  of  the  fluid  in  its  neighborhood  greatly.  Unis, 

there  is  m  interaction  between  the  water  end  the  body  saving  through 
* 

it.  Two  approaches  to  the  problem,  naosly,  the  so-called  "Hjin-Ship" 
theory  end  the  "Slender-Body"  theory  or  "Strip"  theory  fcsw  been  made 
recently  ,  H*  foKaor  was  developed  and  «o@d  by  S&dbell,  Bstera  and 
Stoker,  Hewaan  end  others.  ‘At  employs  a  systematic  partustoaticu 
expansion  in  terms  of  perturbation  fsetoria  in  finding  a  velocity 
potential  while  the  boundary  conditions  sure  respected,  together  with 
the  condition  that  the  bean-length  ratio  be  small.  Zt  la  basically 
a  thres-dimnslonal  problem.  In  the  strip  theory  the  body  is  also 
supposed  to  be  elongated,  but  also  the  draft  is  supposed  to  be  email.. 

In  that  case  the  flow  around  this  body  is  not  considered  three-  but 
only  two-dimensional .  lbs  ship  is  divided  into  vertical  cross- 

For  surface  ships  the  Froude-Krylov  hypothesis  is  satisfactory  since 
buoyancy  is  the  main  exciting  force.  In  the  semi submarine  the  hydro¬ 
dynamic  forces  are,  however,  primary. 
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sectional  strips  which  are  treated  separately  frco  each  otbsr, 
neglecting  tbs  effeato  of  the  loRgitudiaai  perturbation  valoeiti«B. 

The  t cam  and  saaaent  srs  found  on  each  elemnt,  ©ad  on  intagrati  on 
over  the  length  gives  the  total  force  «nd  uxment  on  the  ship.  Die 
advantage a  of  tills  aetbod  over  t fas  tbtn-ahip  theory  ere  lmedletely 
seen.  Qsx*  can  use  the  sisral.es*  two-diSGsasionaO.  potential  theory. 
ltafortun&<5elya  espirleal  correction  f ©stars  have  to  be  ueed  to  rs*S® 
the  theory  correspond  with  eitperisBeatal  results.  The  thin>ehlp  theory 
respects  till  physical  boundary  coMitio ©a  but  has  large  esssputabianai 
disadvantages  and  boa  not  yet  boon  sufficiently  dewslopswa  to  get  a 
complete  picture  eg  the  notions  of  a  Ship  in  yaws.  Saarsforoj,  the 
atrip  it*  ary  was  chosen  to  find  scat  of  taws  hydrodyosa&c  grcpertles. 
Viecooity  of  the  vater  and,  therefore,  friction  e®a  b®  proved  (eaa 
App.  V)  to  bsv©  little  effect  oa  tfca  motioas  with  relatively  low  Frouda 
auebere.  D»y  have  been  igsored  &a  tMs  tsreatEsesst. 

c 

f 

Express!  otto  for  the  heawe  ©setting  force  esd  tine  pitcMsg  exciting 
aeaeat  on  a  submeygsd  body  hem  feeaa  developed  by  SS^lsm  (Hof.  5),  She 
heave  force  in  astern  seas  iss 


5^  in  r*  -t^v-cj^ttpe 


where  a  is  the  wave  eqpiitude,  t  the  free  surface,  and  A(x)  is  the 
transverse  cross  sectional  ansa  at  a  distance  sc  frosa  aidehip.  This  la, 
however,  for  a  body  without  a  strut.  To  ta tea  the  strut  into  account, 
the  expression  has  to  h*  aedified.  TJaa  hydrcdynasle  pleasure  on  the  body 
is  given  by  f annuls  [l4j  of  Raf.  5s  .  . 

jo  *.  fya  <€  ^  |  5W$  V 


where 


and 


P  =»  y-  [*  i-CV-cJtJ 

<5?  ,=-.  -~  ^  c-es  JbL[ic+CV- c)t\ 


Tfce  total  force  in  the  vertical  plane  expressed  in  complex  form  is 


then  given  by; 

*Yl  .42, 


=-  —  ft  j  -f>(c^s6>--  i  si^. 


Therefore ,  we  get:  /5?__ 

—  —  ftj f>  4.$  —  *L<& 

°  ~3r'ty 


which  is  the  heaving  fores  for  a  strip  of  unit  length  at  distance  x 


With/” »  O,  hence  without  the  strut  w(st)  «  0)*  this  give»s 


The  difference  between  the  last  two  expressions)  coes®b  eft® r  integration 
to  about  15  per  cent  less  for  tSss  body  of  revolution  with  a  strut. 


This  expression  is  numerically  evaluated  for  two  wave  lengths  in  App.  Ill 
(pg*  48). 
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A(x)  la  ibis  aggression  twees  the  areas  of  tbs  sector  CAOBC  *t  x. 
Finally,  Zjbeccoee  a  *ua  of  two  terns: 

zl  *  fi  e0*  *  f £  8in  “*• 

Tbis  force  ves  deviated  for  wevalong tbs  of  four  and  olgfet  ft. 
and  gives  u  positive  value  (upward  f caress)  vben  tbs  w&vo  trotigh  passes 
ibe  middi®  of  tfe©  sMp.  Share-fore,  tbs  uafewefcsr  sxsrt  will  counter¬ 
act  ■Mfe®  buoyancy  fares  due  to  tbs  part  tbat  pierce®  tbs  surface. 

The  equation  of  tbs  v&m  profile  is  esmssed  to  bo  eisrausoidals 


y  =  a  ®ia  -£«L  (x  »  et)  ubss©  x  is  an  absolute  abscissa. 
Balstiva  to  tfee  moving  ahlp}  we  tbaa  gets 
y  *  a  sin  ■«-  [x  -  (c  -  V)t  j 

g  Y? 

=  a  sin  — j—  +  (V  -  c)t  j 

She  difference  la  dieplaeacent  la  tbe  save  sjhS  is  csJje  *?at®r  i«j 


[<  +  (V-cJfcJ  ^ 


Sloe*  the  strut  le  eyonetrical,  the  tore  v  sin  -2~~  t*  odd  end  toes 
an  integral  ■*  0.  Besot,  the  buoyancy  force  beeaxutsi 

'T 

The  total  heaving  force  ie  then  z  r  z  u.  *  f,  so*  4  +  rs  tin  arts. 
This  is  also  numerically  evaluated  iirApp.  Ill,  pg.  98.  1 

U»  exciting  aoeoent  due  to  the  waves  consist#  basically  of  four 

p*rtSJ 

1.  The  moment  os  the  body  of  revolution  due  to  tta  orbital  uav® 
velocity  ml  to  the  wave-body  interaction. 

2.  .The  acoeat  induced  by  tbs  change  in  buoyancy  elans  the  length 
of  the  strut  due  to  th»  wav©  slope. 

3.  The  effect  of  the  aww  force  on  the  strut.  This  force  must 
be  tahen  into  account  since  the  point  of  application  of  this 
force  is  relatively  far  fresa  the  origin  of  coordinates.  A 

moment  is  gsasareted. 

4.  The  hydrostatic  hcrinontel  pressure  caa  to®  strut  which  varies 
according  to  the  position  of  the  wave  with  respect  to  the  strut. 

Kaplan  gives  far  the  accent  00  a  near  surface  body  of  revolution 
in  astern  waves  (Ref.  5): 


Although  this  expression  was  develops  tl  for  completely  submerged  bodies, 
^il^m^bs^se^^t^  sufficient  accuracy  in  cur  cose.  Indeed,  dws^effect 
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of  the  strut  oa  the  pitch  nceat  exdtetion  of  the  aeia  tody,  it 

negllgliae  since  the  greeter  pert  of  tfcet  acesrnt  Is  contributed  by  tb 

ertrtmtties  of  the  mein  tody  itself.  lbs  direct  effect  of  the  strut** 

»fcich  lei  celled  Mj,  is  discussed  subsequently. 

L  SevelGp^u  liTods 
£» 

*  ss.  </* 

^£2  4*  s^^(V-clt 

u  A  A  _ 

«s  ^  Ape  Los  f(ujt 

;_U  '  *T" 

_  ^  /V*J  =iK  5I£.  4*  Siu  t 

2-  A 


fe) 


pg.  49 

*  »•  cos  u*  +  h£  Sin  Wt  (see  *es>.  n^.  in  tMs  egression, 

coa  <vfe  is  found  to  be  the  doainsat  term.  Etoace,  we  met  thet, 
since  2,,  the  term  la  sin  <4,,  vne  drasinant,  the  g&nee  angle  between 
beevlog  force  acd  pitching  mceaat  is  sround  9c0. 

2J»  effect  of  changing  buoyancy  can  be  expressed  mt 

sl^  *r[*f  fV-cJb]  d% 


u 

A 


cn  *3T(V-c}fc 


=  -*fj, 

since  the  strut  Is  tywaetrlcal 


L 

T 

L. 

'  4 


2: 


^  =»  a£  eoe  cut. 
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Tb#  third  rnooent,  to  taka  Into  account  is  the  direct  effect  that  the  hyaro- 
dynsodc  pressure  at  the  wave  exerts  on  the  strut.  Hue  pressure  on  fore  and  aft 

parts  at  the  strut  is  not  always  balanced  eo  that  a  hawsoaic  surge  force 
ceueea  a  naeent  on  the  vessel.  ifla  eatuae  the  now  arc uM  the  sstrut  to  be 
two-dieBaeloaa.1  m  thm  Laplace’s  oquasica  as  valid. 


3  y  3  k 


where  ©  is  the  velocity  pofc»s$iai 


F1"- 


®w  wlcseity  potmtiial  of 

&  c~*K\ *  .Ul 


-  at r?  cos  h: f(Y - c) tj 

and  the  hosrisoatal  easqKxssat  a?  the  aeb&t&l  velocity  e.t  depth  h: 

v  -  *,  sr[itf(v-cjt]  . 

Actt^rdiag  to  Bernoulli  ‘a  theorem, 

H4+K=^J 

J2 _ T>i  1  a 

~jp*  —  -g-£  —  j-'V  if  w*  change  tha  velocity  potential 

aueh  that  the  term  in  t  is  iocaarporatsd.  JSL.  msm  the  psrtial 

absolute  St  .  .  , 

derivative  vioh  resspect  to  t  in  anisxls  system,  »£=«*-  v  uw%.  y  ^ ii i 

S5t  dt  da 

*  *  f  +  V  t&ere  £  is  an  absolute  a'  seises.  Sine®  the  strut  is  vary 
aarrov,  we  msy  negxact  tbs  intersctiois  tetass©®  mv®  mA  body,  slsiefc  matyg 


+r'r 


1  - 


if  m  cfearngp  t&J  velocity  potential 


the  potential  of  wave-body  flow  »  o.  Eases ,  &®  taa  fessaula  for  tbs 

pressure,  te  s*sy  put  j)  «  ^ 
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5tHs 


Ml 

A 


[a:  +  (V-cJ  tj 


The  re  mating  surge  force  is  them 

.  .  .  2.  _  zitA 

X,  =  *. 


L- 

2- 


A 


si.  A  f  i>  It* 


4(pc  £- 


'O 


fk 

2. 


“t 


ii 


2-F 

A 


[xf-tV-eJb 

£tr, 


Stud 


Mj  a  X  times  the  Jicaent  arm. 


»-  1 

5l-H  SlIao^  J  CcPS  V"^j  fc 
(see  App.  Ill,  pg.  4-9) 


To  find  the  moment  arm,  the  center  of  pressure  has  to  be  given. 

Since  this  pressure  decreases  exponentially  with  depth,  the  center 
can  easily  be  found  (App.  iffy.  We  get  Hj  »  m'"  cos  -*»-  (V  -  c)t  . 
The  hydrostatic  preswre^!lnSt£ethorizontal  direction  will  have  a 
resultant  different  from  zero  due  to  the  wave  profile.  This 
pressure  (perpendicular  to  the  hull)  vill  be  p  g  t]  and  is  a  function 
of  x.  In  the  direction  of  the  longitudinal  axis  the  pressure  ever 


4 

- 

f 

r4^  il 

— 4- 

dx  will  bei  p  |  p  g  dT|  -~£-  a*  dT|  . 


where  Has  sin 


2n 


x  + 


(V  -  c)t] 


P  *  P  S 
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Hut  total  aurge  force,  when  the  wave  nodal  point  1®  at  tbs  cantor 


otf  tba  body,  la  than? 


\  =*  P  8 


1 


0 


Tp  -Si.  a». 

a* 


nils  is  tba  tin*  of  maximao  surge  faros.  In  this  expression. 


2ns 

T1  *  #ln  — — la  tins  tbi®  fores  trill  change  hsKaaaieelly  ao  that 

A 


we  find  after  a  cat  calculations j 

)C^  ft*  +  z  (  C.<?s  t iu^oC  d%\  cos  ~(\l~cjtr 

'$.)  L  Jo  A  A 

=  c-os  £!T  (V-cj  t 

A 

w  is  the  nsxlouEi  width  of  the  strut.  Xk  is  faud  to  1>»  s  vary 

2SC2C  ®t 

anall  Qusntity  which  can  be  ignored.  ifce  total  jsssoat  acting  an  the 
somlwibr* trims  will  then  be: 


M  ss  (n^  +  n£  +  b^*)  cos  uit  +  ss|  aln  ast  a  ra^  cos  aA  +•  sin  tut.  {See  p.  50.) 

Evaluation  of  Bfrdrodynsslc  Pgrtvstivga 

Kcrvin-Krouiovahy  (Ref.  6  m&  7)  has  allied  tba  atrip  theory 
to  find  the  coefficients  of  tba  aquations  (1). 

Ttm  potential  of  tba  flow  around  a  cylinder  moving  perpendicular 
to  It.  «rt.  i.i  |..f  ...errl.tl  ■*““»  of  «te 
cylinder,,  end  R  Is  a  vector  to  any  point  In  the  flow. 


Cd  tire  surface  at  the  cylinder,  rsB, 
sad.  I)  a  vr  con  cc  .  If  sc eacd- order 
tesm*  are  isaglected,  «  »  o  -li..  and 

3t 


a z 

dx 


T7 

*  2r  j  p  coa  a  dec. 


33»  vertical  velocity  of  the  strip  in  v  a 


sc©  +  V9»  Carrying  out  the 


differentiation  and  integration,  ^jeveitlcsl  foroa  oa  the  'body  induced 

is 

Jsy  the  body's  own  motiorif found  by  ftorvin  to  be; 


=  -  pTT  ('-L  -  yC  9  +  VdJ  k* 


-h 


+  irvL'1  b  OJA/  -  006  4  V  &J  dp 


vhere  tan 


*x  = 


/P 


dr 
d e 


•  Frca  this  equation  follow  5 


r  r~  dx  * 


£1 


h 


ix  =  -  i.  ffrV Jr  tan  ft  d<ic 
d-x  =  -  pTT  |  -t  V  ^ 

=  ^p-rr  V  J  ^  xftwY  Jrzp  dp 

-  1  fnr  Vz  j  f  dx. 

Use  moment  equation  ie  chts3ned  by  mltlplyiag  ewsry  asparasaicn  cinder 
the  integral  sign  in  &  with  -  x. 

4  «+^TT 

Bt  «  -  jLpirV  Jrzpdp  -  -t^rY  j'uia^/b 
C-i.  ==.  +IpirV^  -t-  ^  ^ 

bjL  IT 


P 

Ci  =■  ■+• 


^pirV 


use  subscript  1  designates  the  hydrodynamic  part  of  the  derivatives  a 
equation  1  an  estiaated  by  Ibrvin-Kroukova^y  (Ref.  6  and  T>* 
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In  tills  report  an  attempt  is  made  to  correct  the  foregoing  coefficients 
for  certain  hyorcxiynenic  effects  not  fully  considered  In  JCorvin’e  earlier 
works.  Sheae  include  the  effect  of  frequency  of  notion  and  cf  the  free 
eurfaoe.  In  order  to  distinguish  the  coefficients  so  corrected  from  the  JCorvln 
coefficient#,  the  corrected  coefficients  ere  given  without  the  subscript  1. 

1.  Added  Nms  sad  Added  Moment  of  Inertia,  a  and  A 

It  has  been  found  by  experiment  (Pat.  10  and  11)  and  theory  that  tha 
added  bbss  is  not  a  constant,  as  calculated  by  atrip  theory,  but  depends 
on  forward  velocity,  frequency  of  motion,  sad  also  very  strongly  on 
relative  subtaargenee  for  eubaorgefi  bodies.  At  u>  «  0  the  added  mesa 
theoretically  goes  to  infinity.  Uracil,  ©song  othars,  calculated  the 
change  in  added  mass  from  that  feusd  by  atrip  theory  to  to  both  surface 
effects  and  frequency  of  notion.  &  also  calculated  these  changes  for 
body  shapes  differing  from  a  cylinder.  Zb  found  a  correction  factor  of 
infinity  for  w  ~  o.  Based  on  Ur  sell*  a  work  end  cn  the  results  of  the 
experiments  of  Hof.  10  and  11,  various  values  cf  added  mss  end  maaenta 
of  Inertia  ware  selected.  Ibese  valnee  and  values  of  the  other 
coefficients  of  equation  (l)  sslsated  m  described  in  the  following 
paragraphs  were  than  utilised  in  an  iterative  computer  program  in  an 
effort  to  achieve  maxtiman  correlation  with  the  experiments  of  Ref.  3. 
ra»  beat  correspondence  was  obtained  with  the  added  mass  equal  to  about 
60  per  cent  of  the  actual  mass  at  V  «  0  up  to  200  per  cent  at  the  speed 
where  m  a  0  {V  =  c).  less  drastic  corrections  were  found  for  the 
added  moment  of  inertia.  Hsa  numerical  values  cf  a.  A,  m,  and  J  can 
be  found  in  App.  TV  (pg.  51). 


•  1.6  - 


2.  running  Coefficients  (p  and  B) 

In  a  perfect  fluid,  dissipation  off  energy  (which  i# 
damping)  can  only  occur  by  tbs  process  off  energy  being  carried 
away  by  gravity  waves  at  the  surface  and  through  energy 
radiation  below  the  surface.  Since  the  effect  off  the  free 
surface  was  not  considered  in  the  expressions  of  and  on 
PR*  14,  they  are  not-  valid.  In  an  attempt  to  take  account  off  free 
surface,  Korvia-Eroukovsky  replaces  then  respectively  by  ju(x)  dx  end 

J  H(x}  x20x  +  2VD  -  2Vprr  j  r  tan  p.  >tadjt  where  K(x)  represents 

the  dumping  force  psr  unit  length.  Holstein  (1936)  equated 
H(x)  to  pg  A2p,  where  A  is  the  ratio  of  the  amplitude  of  the 
waves  radiated  by  the  strip  over  the  amplitude  off  heave  off  this 
strip  section.  Ursell  (195*0  ®ad  Grim  (1953)  have  found  values 
that  correlate  to  e  certain  degree  with  exper  lisente.  Hcwaver, 
these  theories  cannot  bs  used  in  this  case  since  they  apply  to 
surface  ships.  - 

A  way  to  find  the  Assaying  would  be  to  consider  the  body  off 
revolution  as  a  hydrofoil  on  which  a  lift  is  generated  due  to 
the  heaving  and  pitching: 

dL  »  -i-n  £2km.  «  v2  AS.  dS  is  the  horircntal  projection 

«rf  a  0 c 

of  a  cross  section  ■  2r  dx. 

<52  »  p  — JL  *  Vs  r  dx.  (|0.) 

d  cC 
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Jet 

Ow  whole  difficulty  lies,  however,  1b  finding  since  It  varies 

over  the  length  of  the  body  for  our  peculiar  shape  with  Its  vary  low- 
aspect  ratio.  Therefore,  this  method  was  also  discarded. 

The  method  selected  in  this  report  far  calculating  the  dsnping 
coefficients  was  given  in  a  recent  paper  by  Osman  (Ref.  8).  Damping 
coefficients  were  computed  for  an  oscillating  ellipsoid  with  forward 
ffiOvica  near  the  free  surface.  The  ellipsoid  ie  represented  by  a 
distribution  of  singularities,  nearly,  steady-state  and  oscillating 


dipoles  and  quadripoles.  The  temping  coefficients  are  then  four/ 
from  the  energy  radiation  at  infinity.  Using  the  sysfcols  of  Ref  (8).. 


the  diaensional  damping  coefficient,  b,  of  tills  report  Is  expressed  ar 
follows:  *1  .^3 

*.«  65J  =  -  BZ /> a.W  z  I  t  X 

40  f  I  +  477  C£>S  it) '/a 


X  (Dj  (u)  /  5^  (-u)  J-U.  {Formula  31,  Ref.  8) 


and 


is  the  nondlaensional  heave  damping  coefficient,  and 
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eare  Green 9 ■  integrals. 


*1'  *2/  *3  ar*»  respectively,  the  half -length,  half-beam, 

and  half -height  of  the  spheroid. 

j{q)  =  spherical  Besaal  Function 
Kn  =?  ellipsoidal  coordinate 


$ :  XW-tn^rf^ .  ^)siAf/2 
t  -at  3 


U  = 

o 


0  <for  T 

for 


Ce*~'  (tf) 


r  > 

similar  expression  was  found  for 


v« 


where  B55  is  nondimensional  pitch  Ramping  Coefficient. 

An  IBM  704  ccKputsr  program  exists  to  find  these  value®  as  a 

function  of  t,  «~»  cad  e^,  ag,  &y  ’*  have  essusad  that  the 

shape  of  our  body  of  revolution  doss  sot  differ  appreciably 

from  a  spheroid,  such  that  the  ssmn  formulae  are  applicable . 

Fig.  2  above  the  speed  and  frequency  dependence  of  b  and  3  for 
different  wave  lengths.  It  io  seen  on  Fig.  2  that  unfortunately 
minims*  damping  is  associated  with  the  resonant  condition.  For 
Minwweieal  values,  see  Ann.  Ilr  (|jg.  52  -  53)* 

3r  Coupled  Banging  Coefficients ,  e  ana  E 

In  an  attempt  to  account  for  the  free  surface,  Korvin- 

Kroukovsfcy  also  codified  the  expressions  of  e  and  E. 


e  »  J  H(x)  *to  +  kgi «j,  i,j.v 
E  *  |  H(x)  x  <5x  -  Icglc^  »V. 

where  kg  and  an  ccorection  factors  due  to  shape  of  the 
body,  end  the  free  surface  (the  average  kgk^  amounts  to  about 
75  per  oent  for  surface  Alps).  Havelock  (R«f.  9)  found 

* 


at  DTMB. 


nearly  the  sane  expression  test  a  floating  half  spheroid: 

e  •  pr/V,  and  E  ■*  -prff  vith  p  =  0*51'?  fear  a  Xsagth  dissetar 
ratio  of  8.  Afekowits  (Ref/0  finds: 

e  =  -(2^  +  VZ^) 

S  =? 

.ter*  the  expression  for  e  is  siailer  in  form  to  that  suggested 
by  Kowin.  Ve  also  see  tint  B  ^  e,  and  this  has  also  been 
proved  by  experiments  carried  out  by  Gerritazoft  (Ref.  13).  e 


was  found  to  have  a  different  sign  and  a  larger  value  than  E. 
To  find  an  expression  for  oca  can  use  Partoula  (10)  of  pg. 


■jiz  - 


£.  a=  _  p  V J  j|&ju  -f  ‘iQ'fm. 

is  fouafi  in  the  cone  ways  -  p  v  j 

'  ~x 

E^perinwats  have  been  conducted  at  ETMB  to  find  for  different 

types  of  deeply  auteerged  outosrinea.  an  following  empirical 
foomula  was  found: 

%  =  0.871  (St,  -  ^Jb'  =  'to**  o'  ~  ~gr  »  ead  Eg  ' 


and  kj.  are  virtual  maas  coefficients  \&ieh  can  be  found  in 
Iarib'a  "Bydrodynamics”. 
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Tfe#  valuao  far  •  and  S  finally  selected  were  based  upon  exactly  the 
asms  procedure  as  that  outlined  on  jg.  3.5  for  determining  a  end  A, 

Shew  ▼slue#  are  given  in  App.  17.1  (pg.  53).  'She  values  e  end  £  found  in 
this  way  are  not  ouch  different  from  those  calculated  using  the  method 
described  in  the  previous  paragraph.  Xu  any  event,  the  values  of  e  end  E 
for  the  bare  hull  body  are  very  staell  coopered  to  the  e.^,  end  B;i  found  for 
the  sendsubmarine  equipped  with  control  fins  st  the  extremities  ef  the 
vessel  (set;  Chapter  W). 

4.  awias  Constanta  c  and  C 

U»e  coefficient,  c,  exists  because  of  the  presence  of  the  surface- 
piercing  strut  and,  C,  because  of  the  metacentric  stability  of  the  semi- 
submarine.  From  equations (1) and (2) it  is  seen  that~c*Z„  and  C»Me+ VM  .  The 
heaving  foroa  due  to  a  aUtergemce  *  is  S  =  -pg  »  «pg  A^i 

therefor®,  e  «  pg  A^«  a pg  A  .x  «p  g  &  *£«■  gl  where  X  f,a  the 
longitudinal  raoesat  of  inertia  of  t&a  waterplaao. 

At  the  M.X.T.  Towing  Tank,  tests  sore  carried  cut  to  determine  the 
hydrodynamic  acraaat.  on  the  aemlcttflawiarim  with  different  speeds  (Esff.  3), 

At  every  speed  it  was  found  that  the  vessel,  once  given  ©  small  initial 
trim  angle,  took  on  a  certain  larger,  stable  trim  mgla  whi oh  was 
variable  with  V  as  shown  in  Fig.  4.  At  the  stable  trim  angle,  the  longi¬ 
tudinal  stability  moment  is  eqjunl  to  the  hydrodynasie  rosssat.  For  smaller 
e  the  hydrodynamic  moment  would  be  greater,  and  for  higjsr  e  the  stability 
moment  would  dominate.  VSills  the  assumption  that  the  stability  moment  is 
epprcocinately  linear  with  8  is  reasonable,  it  is  unreasonable  to  assume 
that  the  bydrodymmic  vxmnt  is  also  linear.  A  parabolic  function 
relating  this  moment  to  8  web,  therefore,  developed.  Shis  is  given  is 
*»•  &  (m-  54)* 
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5*  O*  Ccta&M  Virtual  Wbae  and  jfaaggfc  of  laartia,  6  sod  B 

tese  coefficients  sMcfc  vara  found  by  strip  theory  to  bs  actual  to 

«pj  r%xds  vith  ecrvsotios  facta?  &^9  Om  to  the  surfac©  (see  pg.  l8)> 
are  very  staalX  gaaotities.  as®  easpsriiaBBts  esc  Ref.  lb  also  shoasd  this, 
because  te  instruaeats  ware  mot  mas&t&vz  ©aou^i  to  give  trustworthy 
results.  It  is  possible  that  these  eeapUag  term  also  cheaga  eocsiderebly 
at  very  lew  fragueseies.  Sas  megEitea®  tee®  coefficients*  ©aglecting 
free  mirfee©  and  fsc-guaacy  dstpeaSoney,  eg-©  gluon  in  tlT,  pg.  55* 

6.  She  Couptoa  goriza  Caastxmts,  a  aad  g 

Lite  the  cGOfficieats  c  se&  C-,  t&oss  coefficients  exist  bseatsaa  of  the 
presence  ttf  the  strut.  Froa  equations  (l)  ®a&  (8),  ws  finds  g=  +  V2^ 
aafl  G  a  Kj,.  Sfervi^lQfouitovalsy  (Ref.  7)  tomiZ  the  smm  &K@s®mlam  in  the 

fora:  g^pgjB^x^dx^Vb  G  ^  6  /  B*  k  dx  ■tesa  B*  is  the 

to®m  of  the  strut.  Eacaus©  te  strut  is  spaaetrisal  uitb  respect  to  the 
erosa-eGater  Mao*  Z&  »  0.  &®4rieal  sraletloaohipa  fos?  \  dswloped  at 
OTI®  for  fioep ly  atosrgei  cubeor inss  give  the  fc&lewiog: 

J? 

Zu  «  »{tL  +  S’)  «  L*  a*)  tte  dScxmsiooless  lift. 

w  l/2pLsV 

jj 

De  =,  .... „.,...,..v?.  ia  a  dirassasdoaless  drag.  »  .23b  ffi'  °*^  where 

P/2  lAr2 

2y 

m*  a  **£g*‘  . 

Ojerefcr©*  g  a  l/spv2!,2^  +  D*  )•  $h*®e  will  be  a  sasH  coefficient  G  <  lua 

v  2ha  staple  expression  fcpS^,  is  ussd  her®  in  lieu  of  1 1»  meant 
sophisticated  expression  for  fc  on  ps«  17  **  lBt  because  tfe®  coupled  sprang 
constant  g  is  a  relatively  insignificant  one.  Ibus*  there,  ves  no  nfitrfl 
to  introduce  free  surface  end  frequency  dependency  corrections  hare.. 


to  tha  fact  that  the  strut  does  not  stand  In  the  Biddle  of  the  body. 
Ttm  moment  with  a  submergence  a  will  be;  -  -pg  x3  r<  where 
is  the  horizontal  distant  between  the  center  of  body  and  the 
canter  of  strut. 

\=(4*“)  *  "Pa^s^ 

1 3*  4  =  o 

G  =  +pg  x8  Ay. 


For  the  numerical  values,  see  App.  SVf  pg.  55* 


lit.  solutioh  or  the 


The  equations  of  motion  ftre  of  the  types 
(/Iftt^Z  +  iz  +  CZ  -f-<U9  4&0  2S.  ^C&&  (Jjt  +'^2  SLlitub  (Heave  Eq. ) 

+B6  +Ct»+C.,l  +  J>*  j-Ei+6z.=  'to,  COS  art  -f  TTh  (Pitch  Eg 

Mf 

»  *  (v  -  e) 

A 

The  coefficients  are  either  eonstasat®  as?  functions  of  V,  Siacs  the 
equation  of  pitch  ia  not  line®?,  the  sis^la  mstfeod  of  solving  a  set 
of  linear  differential  equations  cannot  b©  us©S.  B«aves?#  to  maSSse 
a  preliminary  estimate  of  fees?  Use  theoretically  cesaputed  responses 
at  heave  and  pitch  will  correspoad  xrtth  the  experimental  results, 
the  pitch  equation  aae  ma&s  linear.  Xa  this  way  it  is  possible  to 
write  a  computer  program  that  solves  -ass  aquations  ia  a  very  short 
time,  such  that  coeffieiente  that  had  to  bs  found  by  trial  end  error 
method  could  easily'  be  varied. 

The  experimental  pitching  responses  were  used  to  find  values 


for  C6  +  CjT^ejfor  different  sjeeda.  Hhea  these  values  are  divided 
by  e,  a  coefficient  C'  is  feund  which  will  ba  the  coefficient  of  © 
in  the  linearised  pitching  equation*  For  speeds  above  the 
experimental  range*  an  angle  6  was  assumed  and  later  on  corrected 
by  trial  and  error. 

Tbs  solution  of  the  aet  of  too  linear  equations  is  given  by 


Korvin-Kroukovaky  (Ref.  6)  in  the  for mj 


QR  -  FS 


~  SR  -  HP 

©  <*  mmm  saw 

§R  -  P3 


•  2k  ~ 


where  P  ®  -(*  +  *)u£  +  ibtu  +  a 
Q  +  l*us  +  g 

B  +  l&i  +  0 

3:-(J+  A)u£  +  ibw  +  C 

n  ^  rr  ^2.6- <? 
i-i  »*  C  w 

O 

M  »  M  e^M 
o 


.  u 

e  and  <_  sr®  the  phase  angles  between  respectively  exciting  heaving 
force  and  pitching  moment  and  the  wave  motion. 

5  ard  £,  are  tfca  phsw*  engine  betweoa  Mem  saS  pitch  and  t-feair 

(se®  pg.  13) 

respective  origins  of  tints  •  M  «  bs^  coe  wt  *f-  e^>  alls  wt  jh tBSse*  bra , 

6  *  tan”3-  «-SL  and  M  a  /a?  +  s|  .  6  will  be  found  in  th© 

sa  ~  o  V  A  *" 

®1 


farm  §  *  e  +  ito.  So,  ©Q  w./a^  +  b® 


h 

festafi  -£*■.  An  scalagoiio 


remark  goes  for  Z  aufc  A  cotapnter  progress  tses  witlm  tMefe  gives 
the  values  of  -©0  and  4,  and  their  corrasposdiag  pi*®©  lege  for  tl© 
sendaubtaarine  is  fcdlowins  wave®,  4  and  8  ft.  long.  Fig*  3  ®hma 

the  results  of  two  methods  of  computing  the  response  curves  in  comparison  with 
tteexperlaental  heave  and  pitch.  iftuaerieal  results  can  be  found  in  App.  IV. 

yffflSSEP of  solving  the  two  equations  2rlelded  satisfactory  results 
for  the  pitching,  over  the  average.  The  buape  and  hollows  in  the 
experimental  pitch  curves,  however,  are  not  reproduced  in  the 
theoretical  solution.  However,  ta®oc  axprriaental  data  of  both  heave 
and  pitch  for  speeds  of  1*5  ft/sec  cod  up  show  great  fluctuations,  and 
it  cay  be  questioned  whether  they  represent  ths  correct  amplitudes 
entirely. 
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For  the  has;  vice  notions  It  is  found  that  solution  of  the  linearised 

equsfcicaB  *01*5 too  low  amplitudes  for  speeds  teles?  1  ft/sac.  The 
disparity  is  feoffl  50  to  30  pese  amt  in  4-ft.  vavsen  e&A  a  little 
tetter  in  8-ft  t mes  •  Xt  is  possible  that  either  the  struts  the 
surge,  or  a  throe -d3.o*asiGaa.X  factor  haw®  m  Acpoaptaat  influence  on 
heaving.  Thie  iafiuanoe  i»  such  that  it  i ncreacee  tbs*  heave  response 
of  the  seniouteariae.  At  a  of  sfeonfe  1.8  ft/aae,  the  ejtporiaantail 
aad  the  the  critical  data  Join  and  fog*  spsads  theory  gives  toXuss 

larger  than  the  teats.  Oils  is  not  a  serious  deficiency  because  it 
naans  that  any  controls  designed  to  handle  the  theoretically  predicted 
excitations  will  overcontrol  in  reality. 

Vfe  have  no  my, of  eheeSsiau  tm  u©31  tis®  tfesesy  zfapmesntQ  the 

©ctual  motiaia  for  speeds  sfoove  2  S%/mo  (Fx  ™  0.35}.  St  Is  sects  that 

the  lEotiGHG  fnersesae  greatly  and  wall  teyead  the  liasas*  rsags  90  that 

the  liaasarisfid  aquations  So  cot  rspsreeant  tS©  aetusi  jsat&oas  eay  saore. 

The  speeds  of  synchronism  occur  Indeed  at  Freuds  ambers  sterling  at 

0.25  for  tfe©  heave  ia  4-ft  vavao  sss&  eh  0.29  for  the  pitefe  to  4»ffc 

tnn/cg.  These  speeds  go  gradually  up  ssith  the  wave  XecgtSa®.  2a  th£ 

solution  of  the  linearised  equations,  tfcs  paafee  at  synchronism  axe 

apparent;  hut  they  are  only  qualitatively  correct.  Nevertheless, 

the  confutations  show  that,  without  controls,  the  sondsuhnsrinc  cannot 

ha  operated  in  speed  ranges  of  hr  »  0.2  and  up  in  astern  seas. 

linaar 

To  teat  the  validity  of  the /method  of  solving  the  aet  of 
differential  equations,  a  ascend  method,  sanely,  Mara's  numerical 
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method,  vts  used  which  gives  the  amplitude*  of  heave  and  pitch  as  a 
function  of  tiat.  A  cat  of  Initial  cmplitudec  is  given  and  after  a 
certain  transition  period  a  relatively-  steady  feamxaic  action  is 
found.  For  every  speed  the  nsaxlaura  amplitudes  out^toout  wre  taken 
to  represent  the  motions  as  van  also  done  with  the  experimental 
results. 

As  o»a  ba  seen  in  App.  17  na&  from  Fig. 3  ,  tte  heave  response 
in  U.ft.  -saves  follows  very  closely  the  other  solution,  while  ia  8» 
ft.  waves  the  tiae-»hlatospy  ssthod  corraepce&e  mitih  tetter  with  the 

experimental  heave.  The  nonlinear  tern  in  the  pitch  equation  which 
ia  taken  account  of  in  the  T-H  Method  may  account  for  this  improved 
correspondence. 

As  for  the  pitching,  the  cisrvas  stay  close  to  each  other  in  the 
experimental  range  (V  trm  0  to  2  it/ mo)  tilth  the  first  msthod  a 
little  better  than  tba  second  in  4»ft.  waves,  while  it  ia  the 

opposite  case  in  8-ft.  waves  (Fig.  3B). 

For  speeds  higher  than  2  ft/sac  (iV»  0.18),  there  ia  eonoiderahhj 
disparity  between  -tea  two  wys  of  solving  the  equations,  espseially  in 
the  long  waves.  Mara  over,  the  Sisplitadee  go  well  bsyoefl  the  linear 
range  of  the  equations  such  that  they  do  not  corrasporaS  with  reality. 
Although  gusiatitativoly  tfcess  result*  or®  not  valid,  qualitatively 
they  give  an  indication  where  the  high  ssaplitudes  sill  bo  found.  Ohe 
peaks  are,  indsed,  found  at  the  resonance  speed's,  i.e.,  where  the 

tuning  factor  -A  =  «— -  *  1.  For  hmvs  the  resonance  speeds  go  from 

<% 

2.0  to  3.0  ft/sec  for  wav®  lengths  from  k  to  8  ft.  For  pitch  the 


reeooBEie®  fljgoaA  vstsfpi  1®  £rm  3*h  to  k.0  tf'/mc.  It  ie  twn  ia 
fcoth  fee#*®  00$  pitch  mSics®  tiaaf,  tfc&  gsatas  &o  not  coincide  vxmtly 
witli  tbe  fetaratest  psates.  fJMs  my  to  am  to  tb» 

£ mb  ttesst  oSte1  jassaitoses1  tscaa  m&  ti»  sm^p  ssst&om  m y  have  «a 

iaflSSCSEOE, 

Aa  intsseatiag  fesfcur*  of  tte  tite-teisrtOTy  ase^oA  is  that  it 
£&ow3  ©saBctly  Isa?  ®b©  vessel  wmm  in  wssvoe.  St  %jse9  s?o*f  iaetaaaa, 
found  tts&ti  at  mxj  Im  £%®qg&ns&  d?  ©seoantos1  tfe©  wacal  toe  tb# 
teateaey  t-o  oae&U&ts  at  ita  nateal  fjea^usacdes  isiioh  &m  d&ffareat 
fas'  iimm  aaS  pitch. 

Alt-teemed  tte®  SlffostmtlGl  esusv&ieasj,  as  t&a  haw  fomalat®<3  ttessa, 
ero  not  ve2l4  fess*  a  sp»S  rang®  t8s aa  S’s-  «  0.S5*  Utesy  still  asay 

tea  ue®6  if  m  §8©3jsy  mans  to  oosteol  tbs  motion®.  £a  tteis  'Bay  fci» 
eagsllfedaes  oftlS  aitiEsatoly  bs  easaSX  easssgfe  to  safe  tfee  sonatina®  va&i& 
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IV.  TOffi  C0WJB30L  0?  HEAVE  AM?  PITCH 

Ohe  exps-e  salons  for  heaving  and  pitching  motion  derived  in 
the  previous  chapter  form  a  visible  tool  to  Investigate  if  control 
at  these  motions  is  possible  and  how  this  can  be  dene.  We  know  that 
the  equations  do  not  give  the  right  picture  at  speeds  at  and  higher  than 

synchronism  (since  the  response  goat  beyoad  the  linear  s?aoge.  With 
control,  however,  it  say  be  possible  to  de  arenas  the  motions  in 
such  a  way  tlhat  they  reuse!®  in  the  linear  range  bo  that,  we  cea  use 
these  equations  even  for  speeds  close  to  syncteoaisou 

Ibe  aimpXeat  m3  sost  eccmaicsl  way  to  control  s  osrt ©la 
motion  ia  by  ©  hydrofoil  that  can  rot&ta  aloa g  ea  axis  gsrpsadicular 
to  the  streets. 

To  control  tfra  pitching  it  is  edvantagecua  to  pat  the  foil  in 
the  stern  cf  the  ship  since  thia  gives  &  stable  suds,  that  wfesa 
there  is  pitch,  the  lift  works  in  tfc®  direction  opposite  to  tbs 
pitching  angle.  To  control  the  heave  the  moat  favorable  place  for  a 
foil  is  near  the  center  of  gravity  of  the  Ship.  In  this  way  the 
pitching  moment  caused  by  this  foil  will  be  aero.  3aie  prevents  the 
effect  of  the  -nti-pitching  fin  to  be  annulled  by  the  aati-tesevlng 
f in*  Xt  follows  that  in  the  equation  for  pitch  the  effect  of  the 
«atl-plt«hing  fin  alone  Is  taken  into  account,  while  both  foils 
affect  the  equation  of  hs&ve. 

The  wove  is  the  disturbing  force.  The  control  element,  however, 
works  in  this  disticfeeiiee  which,  in  turn,  has  an  effect  on  the  control. 

* 

By  anti -heaving  fin"  it  Is  understood  that  a  pair  of  foils  are  located 
on  either  side  of  the  vessel  in  the  vicinity  of  the  longitudinal  c.g.  of 
the  ship. 
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In  this  my  tee  wri  disturbs  tbo  veswl  «m  tbs  control*  The  rzmt 
efficient  control  would,  la  fact,,  occur  in  calm  water.  By  most 
efficient  it  is  irapllod*  efficient  on  ttm  avsragsj  since,  tilth 
certain  wave  lengths  aad  with  the  veoael  pitching  with  frequency  at 
encounter,  the  orbital  velocity  at.  tea  aafcl-pit«*lr«g  fin  will  bo  so 
as  to  incrssiie  the  effect  of  tea  foil*  At  other  wave  lengths  the 
opposite  would  occur.  Ifeejnsfoa?©,  it  is  better  to  safe©  the  influence 
of  the  orbital  velocity  eras!!,  fflda  can  fey  tea*  fey  placing  tea  fin 


lower  than  the  ship.  She  wave  osfeital  velocity  la  giwsa  by 


and  was  computed  for  different  daptea  and 


wove  lengths.  Sha  change;  in  angle  of  attack  on  ©  foil  dua  to  this 


V 

v_  is  2  -~B-  and  has  feaea  coeluted  for  different  speeds.  It  was 

found  that  etV;  2£t/ssc  for  a  wav©  length  la  8  ft. ,  tee  choaga  in 

*  at  first  judged  to  be 

angle  of  attack  at  the  depth  of  the  csnterXisMs  la  20  ,  which  waaJ^oo 

large.  B»  change  in  angle  of  attach  at  I./lO  below  the  eeatarlino  ia 

10°  which  would  fea  acceptable  f osr  good  functioning  of  tea  aft  foil. 


■Eha  eatl -heave  foil  would  remain  at  the  depth  cf  the  cautarline 


since  th»  pfc&ao  l*g<  between  the  change  of  angle  of  attack,  and 


heaving  motion  is  not  depsndeat.  on  tee  wave  length. 


2ha  angle  of  attack  as  a 

function  of  tins  is  given  by* 

v  L 

CC  zs  ■*"*""  coa  (ulfc  «  IT  ♦  it) 

P  v  1 

coa  (ait  +  it) 

V  1 

-—3k.  coa  (<rfe  +  -y-  tr) 
tdsore  »f4  ia  the  wave  amplitude 
at  depth  h. 


•* 


Depth  of  the  centerline  h  =  .75  ftj;  h/L  ®  .  3.875 . 


3o  " 


To  gat  «a  Idea  about  the  site  of  the  foils  end  the  sanies  they 

have  to  incline  to  control  the  heaving  sad  pitching  motions 

successfully,  one  can  reason  aa  follows. 

The  motions  will  be  completely  annulled  if  the  forces  on  the 

hydrofoils  act  in  such  a  way  that  the  pitching  moment  wad  heaving 

force  which  they  produce  la,  at  all  times,  equal  to  the  exciting 

moments  and  forces  due  to  the  waves,  Ibis  wuM  not  »«  a  control 

have  to 

system  la  the  regular  sense  because  it  voul&jjbnow  from  beforehand 
how  to  react.  In  regular  waves  this  system  could  wor&  but  not  Ik 

a  confused  sea. 

We  shell  first  investigate  the  pitching  motion.  The  lift 
force  on  the  anti-pitching  fin  eera  be  expressed  &s  1/2  p  A  U® 
ubere  A  is  iSfi&S'Sf'ffif'fSL#' WP&P  the  foU' 
produces  is  l/h  A  ifX  for  o c'<  SO0,  c^ia  the  total 


eagle  of  attach.  Jf  6  ~  6Q  ain(offc  »  'tr'j  represents  the  rotation  of 
the  angle  of  the  fin  around  its  am  axis,  than  the  total  angle  of 
attack  becomes: 


The  pitching  moment  will  be: 


■**  -* 

U  is  the  resultant  of  V  and  . 

U*-  +  iYfr  u,s(v,^j 

(V,V  -  mt-X  +  L-jp 
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u3' = y2-  -+  sunlit -f  -UttJ 

Hp  »  Aii:  ^  5tin,(i*>b-z)-  ^cos(a>b^^TrJ  x 

x  [v 1  +  iViTp  $in(u)t  +  ^-Tr} 


His  ubSsnowns  are  SQ,  t,  and  A.  Ttda  mosssat  tees  to  bs  equal  to  the 
©^citing  BESasnt  jb^  'cos  ut  +  tag  sin  at.  xo  simplify  ths  solution* 
the  terns  yd  +  2V  v„  sin  (at  +  ■>«>>  rr)  are  ©scmassd  to  be  negligible 

P  V  A 

coopered  to  U®.  Tills  is  permissible  if  the  level  of  the  foil  is  at 


least  L/lO  telow  the  center  of  toMowsacv  and  fox*  nmeda  hAt&me  tSies 

2  ft/sec.  By  setting  all  of  the  terns  first  in  sic  tut  and  then  in 


cos  at  equal  to  aero  In  tie  equation  Mp  -  cob  at  -  mg  sin  at  «  0  and 
solving  far  *Q  and  t,  it  is  found  that 


(For  a  discussion  on  tbs  magnitudes  of  A  aod  ,  see  App.  VI,  pg.59  'i 

3  « 


6C  sad  t  have  been  calculated  for  an  anti-pitcting  fin  vlth  A  a  556 

SCr 

of  the  longitudinal  cross  sectional  area  sad  -y~>  »  3-4  in  waves  of 


4  and  6  ft.  and  for  speeds  of  2  cad  4  ft/oee.  Itoe  results  (App.  VII, 


jg.  69)  show  that  it  is  impossible  to  control  the  pitching  at  V  *  2 
ft/sec,  while  at  V  *  4  ft/sec  it  is  theoretically  possible  to  eliminate 


the  pitching. 
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ForaiLr.  (l)  gives  the  wxijam  rotation  s t  the  foil  and  (2)  the 

in  order  to  control  the  pitch. 

phase  angle  betssea  the  Motion  End  the  origin  of  tisao^  To  control 

the  hsava,  a©  effect  of  bot&  foils  has  to  be  tafssn  into  account. 

The  anti-heaving  foil,  vMefe  is  situated  midaMgs  at  the  hei^it  of 

the  center  of  buoyancy,  ie  aseunsafi  to  hsv®  the  asm  aMps  esd  are®, 

m  th©  ent:  l -pitching  fin.  So,  the  total  heave  exciting  farce  from 
both  the  pitch  and  heave  foila  la; 

F  a  (|^)V“  4  S^wt-TiJ-  jt+  4lTt) 


v  .  particle  velocity  at 


v  =  particle  velocity  at  depth  h 


6^  —  maximum  rotation  of  center  foil 


n  phase  lag  between  center  foil  md  origin  of  time, 
liis  force  ie  pat  equal  to  the  exciting  force: 

F  =  t\  cos  tot  +  t0  &in  out. 

Putting  the  terms  in  cos  «<t  and  sin  <»t  on  both  sides  of  this  equation 
2 quel  to  esch  other,  it  is  found  after  some  calculations: 
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£05  f.  =  f*  - 

and,  after  replacing  $0  eafi  v  by  their  values  in  (l)  and  (2)(pg.  31); 


ft)  L0$J<a*~ 


*** 


It  is  found  that  for  speed*  of  2  ft/sec,  it  is  impossible  to  control 

tbs  heave  ideally  (App.  VII,  pg.  64).  This  is  similar  to  the  conclusion  for 

-  witch . 

These  computations,  whioh^are  in  themselves  only  apprcalnatlcns  • 
of  the  real  physical  conditions,  chow  that  It  is  Beosatsary  to  take 
the  greatest  area  for  the  foils  possible  and  that  there  will  always 
bo  seam  pitching  sr»d  heaving  at  low  speeds  {Fs.  <J.25). 

A  snore  realistic  method  of  control  is  to  use  the  actions  of  the 
vessel  as  input  to  the  control  elaoaat.  The  control  fins  will  work 
In  such  a  way  as  always  to  counteract  the  motion,  the  vessel  performs 


at  that  moment 


*3^  " 

!U»»  motior  equation*  ««n  then  be  expressed  as  follows} 

^plt  0^2.  ■{■io'i-  -b  (LZ.  B  i  — ^j+Z^^-fa( --(SM "f~21 

Cr*$  +  B4+C0  +  &ztEi+Gz-C..£  • 

~  (&+<* -e-)  +  Mj (t-trk - £))  +Mg( $ f  A 

where  Zg  and  Mg  are  hydrodynamic  derivatives  of  respectively  heave  force 
and  pitch  moment  due  to  angle  of  attack  on  the  foils. 

$  s  angle  between  foil  and  loagltiidisM>X  ea is  of  th®  ship 

sc  =  tangle  of  attaal:  induced  by  the  waves  with  suspect  to  the  horlsontel* 


•  tl 

X  and  o L  are  negligible  quantities  for  waves  as  long  as  or  longer 

k  «V 

thaa  tlKi  ship*  6  am  ft  core  essaXl  fear  tte  vqb&qX  'In  mt®xn  mm  siace 

I  u  *  ’«* 

the  notions  have  a  small  frequency*  Therefore,  ois  ois  <?,  end  6  are 
considered  to  be  of  second  order  and  ere  neglected. 

In  the  equations  of  motion,  the  coefficients  a,  b,  c,  ....  do 
not  have  the  same  value  aa  la  Up  unK^hrolXed  vessel. 

The  added  mass  of  two  foils  placed  at  xn  and  x„  from  the  center 

ir  n 

of  gravity  of  the  ship  is  given  by  (Ref.  l6): 


•ye 


<Lf  -  ~lpT<‘  {  ~L— 

1  1-  V  ,fT’ 


..  <  ** y 

where  a^,  bg  are  respectively  span  and  chord  of  the  anti-heaving  foil. 
Analogously  for  ap  and  bp. 


Zqf  <Hxh  +  **p)  “  “d* 

=+^pTf  (Ex0  +  Pxp)  a  -Df 

V  ®  (EXH2  +  ^p2->=  "Af 

In  our  case  Xg  ■=  0. 


These  veluet.  *>r«  added  to  the  respective  coefficient*  of  the  vessel 
without  fins  (%p„  XV.). 

h*.,  B,?J  e*,,  and  Ef  are  caused  by  the  lift  and  drag  of  the  foils 
due  to  tbs  ship’s  own  raotlon. 

1  -  i-L  *  1/20A 
3  a  1 

D  =  l/S^ACjj 

Frcss  tfess®  two  expressions  ths  teeing  coefficients  can  he  found: 


“_bf 


=  +  1/gpA  V 


l/hfA  V  L 


1/8  p  A  L 


«* 


lA^A  V 


J-fL 

t  S  of 


Computation  of  g. : 

0 

V  a  (Zg)p  6p  +•  (Z&)H  &H 

dp  =  angl®  between  longitudinal  axis  of  ship  sad  the  pitching  foil. 

6g  as  angle  between  longitudinal  axis  of  ship  and  the  anti-heaving  foil. 
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(Zt  )p  &  (Z5)jj  when  both  foils  art  equal  to  each  other. 


(Z*)p  .  l/2pA  V2^-™  +  C„jn  s4  . 


Analogously, 


,2  .  /  ^ 


H,  a  x/4p  A  V*  L  j  -~~2~  +  Cp  j  »  z6  *  1/2 


The  control  angles  and  Sp  sre  mM  dependent  on  6  end  2  such  that 
relatively  simple  controlling  devices  cam  to®  usad. 


8p  »  kj0  +  lye 
*  0 

®H  *  ^3®  +  k42!‘ 

vhere  k^,  kg*  &3>  !%  are  four  control  parameters  to  be  determined. 

If  time  ls$s  in  tha  response  of  the  control  oystem  ara  token  into 

account,  the  angles  becaae  r 

£P  ~  A±0  -  -j-  -i2t4z 

^  *  i5<9  -  tte  +  A^-z  ~  b&  1 

there  t*  «re  the  time  legs  in  the  control  system .  Hence  the  beave 
equation  becomes: 

Jjw'  d.  +  Z,  £$*  +  4.^  )J  5:  4-  [■i+'if  -  Zi  (-fz+  A)\ Z  •+  ^'Z. 

•+  [i  +  <if+2^^i(4JL+4J)Jl9  +  [e  +  «(  -  7$(l.t+£s)]()  +l*f -+2Z^]  0  _ 

-  ^  co s  t  +  f x  n>t  +  'Z.g  I*# h) 


~  —  -  cos 

r  V 


t  L  V 

|uab  +■  —*»»  if y 


is  the  angle  induced  on  the  anti- 


pitching  -foil  by  the  vave.  Analogously,  for  the  anti-heaving  fin, 

-  2£L.  cos  uA.  rdie  total  angle  of  attack  on  the  anti-pitching 
V 


fin  is  then  <*p  +  6p  -  ©• 


When  it  is  fixed  in  the  horizontal  direction. 


-  37  - 


The  pitch  equation  becooetf  s 

[j+A  +  Af  <-[6<6t  -faille  *[<!.+ Ms]9  +Ct  JL 

-e[t+i{  +  A(5fc441]z  +  [E).%_/‘|pf1Ji  +  Q-z.  _ 

—  4n-  t  Co$  U)  t  -j-  cos  to  t  fif  g  c(£ 

The  different  coefficients  of  these  equations  ar®  evaluated  in 

4 

App,  vm.  As  in  ths  uncontrolled  motion,  two  methods  of  solving 

the  equations  of  Motions  ware  used;  one,  the  liaesorixed  method  describe*, 
on  pg.  2A  end  the  other,  Meet:'  lutericsl  Method  (pg,  2 5  -  27)* 

In  l»~ft  waves  (X/L  a  1),  it  is  seen  {A*®.  VIII)  that  the 
motions  are  fsiariy  wall  controlled  if  the  time  lags  are  token  small 
enough*  Large  time  lags  for  high  speeds  can  have  such  aa  effect  as 
to  drastically  diminish  the  virtual  mesas  and  mc*a@Bt  of  inertia  of 
the  chip.  This  effect  is  ©von  stronger  in  the  motions  in  8-ft  waves. 

Thus,  heaving  and  pitching  esa  be  greatly  increased.  Even  an  increase 
in  foil  area  cannot  mates  up  for  the  lose  la  virtual  inertia  (Tables  i'l  &  Si). 

The  solution  of  the  linearised  equations  gives  smaller  fimplitudess 
for  heave  and  pitch  then  Adana’  JteraicBl  J-fetfacd.  This  wag  to  ba  expected. 

Just  as  in  the  uncontrolled  motion,  slnct:  Lite  uonliaaar  term  increases 
the  pitching  {App.  IV).  Therefore,  end  also  because' of  lack  of  computer 
time,  this  method  was  discontinued. 

(Tables  IV  and  V) 

■  Since  the  motions  in  8*ft  waves  were  harder  to  control^  further 
calculations 

were  confined  to  these  waves.  Finally,  the  control  in  4-ft 
waves  was  investigated  with  the  system  that  gave  optimum  control  in 
8-ft  waves  (see  pg.  kl). 

The  numerical  results  and  details  axe  given  in  App.  VI, II,  while 
a  description  of  different  control  systems  that  were  investigated  > 

* 

A  convenient  summary  of  all  the  data  of  App.  VIII  is  given  in  tabular 
form  on  pg.  75a. 


follows  below. 


4  a 

The  system  «  k^e  +  kg«  -  fp  (time  lag) 

*  4 

6g  ~  Js^S  +  k^z.  -  %  (time  lag) 

aid  cot  give  satisfactory  reunite.  Till#  cen  be  explained  not  only 

by  the  tin*  lag  but  also  by  the  fact  that  fear  speeds  higher  thee  •  / 

3  ft/sec  (Fr  «  .26),  the  pheea  ssgbs  between  heave  end  pitch 

becomes  1B0°,  *&ich  means  that  tbs  aati-hsavlng  foil  is  counteracted 

by  the  foil  aft.  To  eliminate  this  effect  the  eati-pltcMag  control 

has  to  be  auefe  that  only  a  pitching  motasnfc  is  produced.  The  solution 

is:  two  anti-pitching  fins,  one  fore  and  one  aft,  which  rotate  in 

opposite  angles  sad  whose  resultant  response  does  not  produce  way 

(App.  VIII,  Table  VI) 

vertical  force  but  only  a  pitching  moment.  From  the  results  can 

be  oeen  that  the  anti-heaving  fin  ia  now  more  effective:  but  the 

pitching  control  is  worse,  since  the  presence  of  a  fin  in  the  forward 

end  of  the  vessel  has  a  negative  effect  on  theSrSfekng  stability. 

This  new  system  requires,  of  course,  alteration  of  many  testas  in  the 

motion  equations  (App.  VXXX,  pg-  70). 

To  sabe  up  for  the  loss  of  stability  In  pitching,  a  different 

set  of  foils  we#  introduced,  namely,  flap  fins.  They  have  a  much 

higher  than  the  conventional  simple  fin.  The  Sperry  Corporation 

manufactures  anti -rolling  gyro  fins  which  have  a  of  1.5  at  an  18° 

angle  of  tilt,  with  the  flap  at  an  additional  angle  of  20°.  The 
d<% 

—■■■».  with  this  type  of  fin  amounts  to  A  further  increase  in 

d  a 

control  ability  can  be  obtained  by  a  "stepcontrol, "  which  means  that 
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the  aegis  of  the  foil  moves  stepwise  free  a  mssteum  angle  In  cote 

oense  to  a  raasiraura  angle  in  -She  opposite  saana s 

v  -  *i  jt 


M 


Is  this  way  the  anti-pitohlag  fall  switefcaa  av®r  wfesn  the  pitching 

motion  readies  its  saalnms  eeplibude  ®a&  reverses  ita  @©z@s.  The 

antl-baswisg  foil  works  analogously.  The  result  with  these  ©arfcmae 
3Cr  , 

cosditioiss  of  -j~=~  and  control  sssthofi  was  an  almost  couplet® 
elimination  of  motion.  In  feet,  the  vassal  is  ,,OT®s?eQafcroll*&>  ** 
since,  the  motions  degenerated  into  aiasosb  a  vibration.  Ks© 
accelerations  showed,  indeed,  very  hi^h  values,  while  ths  frequency' 
of  the  motions  ueo  ccmsidsrshly  Mghsr  than  w  (Table  VII,  pg.  72). 

The  earns  atepcoutrol  with  the  aass  foil  characteristics  was 
used  in  th©  next  tryout  but  with  only  one  s®ti*gitahiag  fin  (aft)  (Table  VIII). 
The  result  stowed  a  fairly  goad  pitch  control,  but  again  the  heaving 
feaeccBs  too  large  at  speeds  of  reso&ssac©  (Vs  3  ft/sec).  Once  ths 


vessel  is  over  t&ie  speed,  the  control  works  efficiently.  However,  for 

these  high  speed*,  the  value  a*  for  «**«&•  Is  rather  high,  sine© 

S  a. 


cavitation  my  beccm  a  problem.  It  is,  therefor©,  safes*  to  take  a 


lower  awism  lift  coefficient  at  those  high  speeds,  at  around  1.2, 
dCl1  .  v^r  . 

which  gives  a  — «!—  of  approximately  k ,  Furtfcosracjre,  it  is  realised 

a  CC 


that  the  large  strut  needed  to  support  "the  large  pitch  foil  at  the  very. 

stern  of  the  semi submarine 


See  App,  IV. 


sty  create  considerable  structural  problem.  Moreover,,  the 
Increase  In  dreg  due  to  thi*  utrut  might  be  so  large  at  hi# 

•peeds  that  toy  resistance  advantage  of  the  aemlsubuarine  over 
the  conventional  thip  could  be  reverted.  Jt  la  alao  true  that 
the  higher  the  forward  speed,  the  lest  the  disturbing  Influence 
of  the  waves  on  the  angle  of  attack. 

These  considerations  led  to  the  introduction  .of  tiro  aatl- 

each 

pitching  fins  of  snail  size  (jL.8  per  cent  of  the  longitudinal  cross 
section  area)  that  era  directly  attached  to  the  for®  and  Aft  aads 


of  ths  body  of  revolution  at  the  level  of  the  apis.  The 
anti-heaving  foil  h m  s  2.J  per  coat  area  end  raa&lae  sziidship.  A 

6Cr  6  Ct 

of  1»  and,  atepcontrol  ware  aasuaisd.  The  value  cf  for  — «-k- 
6  ec  ®  ct 

is  more  realistic,  since  the  lift  coefficient  cannot  be  token  as  lerge 

at  high  speed  a  so  as  to  avoid  cevltation  (sss  App.  VI).  Tm  results  (Table  IX) 

with  this  system  were  aa  good  as  tilth  the  previous  one,  while  the 

laaKliaum  angle  cf  attack  on  the  foils  is  asAllar,  which  QiainlfSiea  the 

danger  of  stalling.  However,  again  ths  pitching  araplitudae  have  a 

tendency  to  oscillate  and,  this  tlas,  not  'wound  e  =  0,  but  around  a 


negative  value. .  Thie  causes  the  vessel  to  m&mB  large  pitting 
angles  (10°  and  more),  which  is  not  allowable,  ©is  saw  was  true,  in 

leaser  degree,  for  the  heave  notion.  Fitch  angle  and  heave  displacement 
signals  were,  therefore,  added  to  the  next  tryout: 

*P  3  *1  J§j  +  *2® 


(App.  VHI,  Table  X,  pg.  73) 


.  Ill . 

It  was  found  that  this  system  Is  vary  effective  to  Unit  the  pitch 

•t  ell  speed*  ,  vfcila  the  beavirg  enplitufiee  reran, In  significant  at 

the  speed  of  resonance.  A  better  result  at  this  speed  is  obtained 

when  the  Jseave  displacement  signal  is  taissn  away  (Ag$>.  VIII,  Table  XI ). 

hfeay  more  systems  of  control  have  been  tried,  but  it  seems  that 

* 

© 

the  system  «_  =  feu  +  iu© 

*  A  /©/ 


gives  the  best  results.  Tt,o  pitch  remains  telow  5°  around  the  speed 
of  pitch  resonance,  while  for  higher  speeds  it  is  much  lass.  The 
heave  takes  values  of  2.5  tines  the  wsw  height  at  the  spaed  of  heave 
resonance,  but  outside  this  range?  it  boeoEss  negligible  (see  Pig.  5). 

Ibis  system  applied  to  the  vessel  ia  Ut  waves  ( X  =  L)  Shoved 
almost  a  complete  dioappaapance  of  any  motions,  even  with  cystoma 
that  still  showed  large  amplitudes  in  8»ft  wavss. 

It  seem*  that  the  control  of  this  ssmiaubsarisas  ia  the  most 
difficult  in  long  waves  which  have  high  energy  an d  cause  frequency 
of  encounter  around  tfes  natural,  frsqusnciee  of  the  vessel. 

Stapcontrol  is,  however,  an  idealised  assumption}  moreover,  it 
is  not  even  desirable  since  it  causes  vibrations  aad  high  accelerations. 

From  the  computational  results,  it  can  be  seen  that  a  harmonic 
control  function  is  sufficient  at  hi#  speeds  (Sr  >  0.35)  l»t  that  step* 
control  or,  rather,  a  "rounded"  atepcontrol  is  necessary  at  speeds  of 
resonance.  It  is  advisable  that  time  lags  in  the  control  system  be  as 
small  as  possible. 


v,  cottauarcwa  m>  RroatsitD maw 


Previous  reports  of  experlinantal  work  emphasised  tbs  small 
motions  of  the  sssisubtasrios  In  ifcead  ness.  However,  with  no 
controls,  ewe  motions  were  eacmurtered  in  regular  astern  seas 
when  jrescmmat  wn  components  yore  encountered.  Sharp  peaks  in 
the  motion  wspiitudc  at  resonance  have  also  been  found  to  ba  true 
in  this  study,  which  is  entirely  sasalytieel  in  nature.  Different 
existing  theories  that  apply  to  either  surface  ships  or  deeply 
submerged  euteasriaes  have  be©a  adapted  and  eg piled  to  the  problem 
of  the  eesi?, submarine  near  the  free  susfffioa.  It  has  bscn  found  that 
the  motions  go  far  beyond  the  linear  roege,  and  the  accelerations 
for  shove  the  limit  of  comfort. 

This  report  shows,  however,  that  with  activated  controls 
incorporating  medium-aimed  control  surfaces,  it  is  possible  to 
reduce  the  motions  to  very  nominal  values.  Since  the  frequency  of 
encounter,  where  severe  motions  occur,  are  relatively  small,  the 
control  machinery  would  not  hsve  to  be  of  a  high-speed  typej 
neither  need  the  poww  requirenect  be  ierga,  since  tbs  hydrodynamic 
exciting  forces  and  moments  are  not  large. 

The  feasibility  cf  an  effective  notion  control  coupled  with 
the  fact  that  the  seoistabearins  cannot  experience  slamming  tends  to 
confirm  the  expectation  that  this  type  of  vessel  should  be  able  to 
sustain  its  maximum  speed  in  very  heavy  aoas. 


To  err iv»  at  a  more  accurate  prediction  of  the  motions  in 
waves,  mare  research  la,  however,  accessary,  Sopeciully,  the 
problem  of  the  added  mess  end  added  ament  of  Inertia  ere 
properties  .that  rtlculd  be  investigated  more  specifically  at  the 
lew  frequencies  of  encounter  that  occur  In  astern  wave*.  In  this 
report  the  values  for  these  coefficients  ware  determined  more  or 
less  by  a  trial-aad-«r.ror  procedure,  which  ma,  nevertheless, 
consistent  with  eoaa  existing  theories  sad  exparissntal  results. 

Tbs  same  remarks  go  fen?  the  coupled  deling  coefficients  whore  the 
existing  theories  ere  not  found,  to  eonrospand  with  experiments. 

Tbs  influenoo  of  surge  and  cures  forces  boa  not  been  studied 
fully  in  this  vest.  Tbs  problem  would,  indeed,  beoaa*  much  nor* 
complex  (see  App.  Ill)  and  forms  a  topic  of  research,  by  itself. 

The  surge  motions  may  have;  a  measurable  effect  cm  the  pitch 
and  heave  of  the  eoaisuboBrine,  sinco  the  large  strut  causes 
coneiderSbbo  moaanta  on  the  vessel.  It  may  be  useful  to  investigate 
the  coupling  of  the  three  motions  (heave,  pitch,  end  surge)  in  the 
future.  All  of  these  reccesaea&atlons  would  mainly  refine  the 
eolations  so  that  the  correlation  with  wxperiaental  results  would 
be  more  cooplete.  They  should  not  alter  the  conclusions  previously 
dream. 

A  natural  consequence  to  this  report  would  be  the  application 
of  methods  and  results  of  this  study  to  the  motions  in  sn  irregular 


gaa.  If  tfca  principle  of  euptapswitioc  la  mtn'mA  u>  be  true,  tlUiii 
AouM  not  poee  way  particular  difficult.!®*. 


APj>,m>IX  I.  CPF3BT3  fl£  3EJg3UBMttOHB 


Body  of  Revolution 


Station 

0 

O 

0 

1 

0.455, 

0.207 

2 

0.702 

0.493 

3 

0.872 

0.762 

4 

0.971 

0.942 

5 

1.000 

1.000 

6 

0.976 

0.952 

7 

0.910 

C.828 

8 

0.778 

0.605 

9 

0.530 

0.281 

10 

0 

0 

Length  st  4  £ t 

Strut 


Station 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


®/®Z8J!8S 

0 

0.253 

0.523 

0.769 

0.939 

1.000 

0.939 

0.769 

0.523 

0.253 

0 


Ijm  <■  £«4  ia 

Vaterplon*  area  A^s  0.24  s$  ft 
Diaplaceonat  of  body  of  revolution  *  0.48  ft 
Piapiaceaent  of  sail  a  0.124  ft3 

natural  heaving  period  a  2.3  sec;,  natural  heaving  frequency  s  2.73  cycles/ 

sec 

Haturatl  pitching  period  -  2.42  sec;  natural  pitching  frequency  =  I.83  "  " 
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Apssmm  n.  gas  imms as  op  smss  aw  arai  mumcm  c&  wncss 

While  the  effect  of  the  surge  exciting  force  on  the  pitch  exciting 
moment  bee  teen  pert  1  ally  taken  into  account  In  this  report  (see  pg.  Il¬ 
ls),  for  the  most  part,  it  hm  not  been  possible  to  introduce  other 
effects  of  surge. 

If  we  call  u  the  surge  velocity,  ve  have  that  the  time  dependent 
velocity  Vo  VQ  +  u  where  u  is  a  harmonic  function  odf  t  in  regular  waves 

and  VQ  is  the  mean  velocity  of  the  ship  in  waves. 


xne  rrequency  or  encounter 

u>(t)  =  “T™  (v0  +  y  "  e 


Mat 

n  tuv/rTj.  VC« o 


®is  experiments  (Bat-  3)  haw©  shorn  that  t-Ma  swgs  u  «  teecsss 
consideratols,  up  to  30  p®r  osat  of  tte  forward  velocity.  2n  etoad 
waves  where  «Ct)  **  (V  +  e),  this  surge  would  not  have  a  large 
influence  but  la  satera  waves  the  eaiotaat  V  »  e  can  teecsss  <pii© 
small  so  that  the  frequency  of  encounter  csa  chesses  ccnsiderebly 
during  the  tints  thet  a  wavs  length  passes  to  chip. 

The  coupled  equations  of  nation  would  becomes 
(a  -  jQu  -  \vl  -  »  Xyg  -  Xg2  -  JCq©  -  (\  ♦  VlQe  -  ¥%e  a  X0cos(wt  -  £x) 


-2£u  -  Z^tt  +  (in  -  2£)s  -  Zyi  -  2.aa  -  2^1  -  (Z^  ♦  VZj)4  •  Z0cos(wt  -O 

<u  -  \n  -  S£s  -  tty  -  Kae  4  {2y  -  1^)3  -  +  V*Qe  -(%  *  Vi^e  = 

l^cosCuA  -  €b) 

X^,  Xq%  Zy”,  are  small  quantities  and  csra  be  neglected.  They  are 
either  coupled  ta©s©a*  or  coupled  moments  of  inertia.  and  Xw  ere 
due  to  added  resistance  from  waves,  dsgesdeat  on  V  «  V0  +  u.  \  is 


increased  drag  due  to  change  in  draft,  dependent  on  Vs .  ss.  £vq  +  u)2, 
\  win  be  the  change  in  lift  force  dua  to  change  in  forward  velocity, 

alao  a  function  of  Vs  -*  (vQ  +  u)£.  From  Chapter  m  it  follova  that 
the  coefficients  2£,  Z^#  !{,  and  end  are  functions  of 

V  «  VQ  +  «. 

Solutions  for  the  three  equations  of  motion  given  an  pg.  h6  can 
only  he  found  by  approximate  methods  with  high-apeed  coeagutera.  The 
problem  is  essentially  a  noaU.neer  vibration  problem  whsra  many  of 
the  coefficient®  ere  difficult  to  define  K&lhosaatieally.  Store  research 
on  this  partieolarsub^sct  in  highly  desirable  to  this  case  of  saai~ 
submarines  but  has  not  been  undertaken  h©r©» 
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APPB'flDlX  m.  KXCrCIRG  FORCES  ASP  HOMEfTIS 
Two  wave  lengths  war©  considered  to  get  an  ids®  of  the  notion 
of  the  eemiaateariaa  la  w®s©e.  They  are  respectively  egoel  to  the 
length  and  tola®  the  length  of  the  vessel.  They  ware  qbcKvm  bseanae 
tbs  natural  fragtescies  of  th®  vessel  «pe  low  sudb  tfest  synchronism 
with  tbe  encounter  frequency  is  likely  to  occur  at  operating  fascia 
for  these  wme  lengths. 

,PS-  7 

The  integrals  in  formula  (UJ^war©  solved  by  Sirapaoa's  Suits. 

Ten  stations  were  considered  sufficient  since  tbe  body  is  abaxder, 
and  the  dimensions  do  not  vary  such  from  station  to  station.  ®te 
results  of  the  coEputstices  for  X  «  4  ft: 

Zj_  s  (0.078  -  0.0008  V)cos  1.57(7  -  4.$4)t  -  (0.69  + 

0.05  V)sin  l.57(v  *>  4.54 )t. 

X  »  8  fta 

Zj_=.  (~  2.813  +  0.15  v)sin  o.t85(v  -  6.4)t  +  (0.13  ** 

0.01  V)co®  0.t85(V  -  6.4)b. 
pg-  9 

Formula  (5^® Ives  for  X  -  4  ft: 

Zg=  +  0.813  aia  1-57  (V  -  4.$4)t 

X  =  8  ft! 

Zg  =  1.17  ain  0.785  (V  -  6.4)t 
Therefore >  the  total  heaving  forces  'become  for  Xs  4  ft: 

Z  =  (0.13  +  0.04  v)sin  1.57(V  -  4.54 )t  +  (0.072  - 
0.0003  V)cob  1.57(7  -  4,j4)t  (ibs) 

X  -  8  ft: 

Z  -  (-1.1  +  0.2  V)sia  0.785(7  -  6„4}t  +  (0.13  » 

0.01  V’)co8  0.785(7  »  6.4)t.  (lbsj 
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The  expression  a  usflsr  the  integral  alga  in  formula  (6)^  were  calculated 
1w  SlwTMort’a  Buies 
1=4  ft : 

(-1.0  +  0.071  V)coa  1,57(7  -  4.54 )t  +  (»o.o4e  + 

0.0003)ain  1.57(7  -  4.54)t 

l  ~  8  ft: 

«(l.605  +  O.Ip  V)eoa  0.785(7  >»  6.4)t  +  (0.(394  *» 

0.020  7)sia  0.785(7  -  6.4)fc 
tepreasioa  ^J^eoaias  for  \~kiti 

Mg=  -0.246  ao3  1.57  (V  »  k.5h)t 

X  ~  8  ft3 

14,  a  -0.261  eoa  0.785  (7  -  6.4}t 

pg.  1 2 

SKpp3B8ion  (slagtve®  for  >.  “  4  ft: 

X3  0.26  cos  1.54  (7  -  4,54)$ 

For  X  =  0  ft? 

=  0.762  QOS  0.785  (V  »  6.4)t 

To  find  tfe©  paint  of  application  of  t&lp  surge  force,  we  assume  that 

the  bcrlxcaatal  hypodynamia  pressure  decreases  expcxisatiaily  with 

2nh 

dapth  according  to  the  law  e”  *■  h.  The  esater  of  preeaura  will 
then  be  given  by: 


For  X=  4  ft?  E  =  0,226  {t 
Foe  X  -  8  ft:  H  *  0.24  (fc 
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5Sus  moment  tfcsres  beccma  toff  1  “  4  ft: 

M3  =  0.26(0.75  -  0.526)cob  I.57  (V  «  4.54)t 
0.136  gob  1.57  (V  »  6.4}t  ft  lbs 

1=8  ft: 

H  ~  O.388  cos  0.785  (V  -  6.4)t  ft  lbs 
3 


pg.  13  ** 

$}  jgivaas  ~  -v  o.O:J  sin  «•»»*» 

a  fare®  that  caa  b®  negLecteft  ,  compared. 

®»  aa©B  remm'k  goes  for  3£^  at  1“8  ft. 


M 


% 


+  M 


(V  ®  4.g4)t  at  k  tMcSi  is 
to  tfes  tQ^Srofi^aesnle  sarge  force- 
Iso  total  moment  then  beeesaes ; 


X-kfti 

M  ■=  {»  0.9  +  0.07  V)cos  1.57(V  -  4.54)t  +  (0.0042  + 
0.0008  V)sia  1.57CV  «*  4»§4)t  ft  lbs 

1  ~  8  ft: 

M«»Cl.538  ■«•  0.15  V)eos  0.785C V  -  6.4)t  +  (0.094  » 
0.02  ¥}i3is  0.785( V  «  6.4)t  ft  lbs 


ft  soft  A  -  Addad.  mm  eal  A’icwi  Unseat  of  Xtwrtifc  (see  $$.  15 
of  text) 

Xa  Ref.  3  tfes  geoantxlcftl  c!*t*.>.  of  ifc®  fwslsatoar&s#  «r©  found  t 


Pbody  *  0.i»8 
P- strut  *  0,07  ft’ 
ra  a  0. 55  is  l«g©  *  1.1  elnga 


J  a 


2y  »  pray8  I/'  Sy*  la*  ®@o&a*  gi'sen  by  8©f.  3: 


fe2  »  0.935 

’■Sherafopa,,  3  «  0.6. 


1, 

“(T 


--  vi 


ss  JL  & 


*  a  « 

«»■** 

p 


H,  a  k  PnJ  i?2 

E  1  P 


2 


ex  =  Is  I 

P  V 


By  aa  iterratlv*  grcxietec  that  e©u<;$at  to  ogtitslto  the  correlation 
vitli  the  eaq^nlasatal  recults  of  Ref.  3s  ^  aod  were  found  to 
bswe  value®  ranging  froa  0.60  at  high  frequencies  to  2,0  at  very 
esiftll  fr*auenelu«  of  encounter.  a»  following  tables  give  the 
final  values  fa?  ia  i»ft  (c»  4, >5  ft/a ec). 


kl 


a  +  e 


+  J 


0.1 

6.97 

1.7 

1.0 

0.7 

6.04 

1.7 

1.0 

13 

4.«)4 

1.7 

1.0 

1.9 

4.15 

1.8 

1.0 

2,5 

3.80 

2.0 

X*  X 

3-1 

2,28 

1.2 

u 

b'M 

iJT.4  wt< 

3.0 

3*0 

i:§ 

4  ..8 

0.41 

,1.8 

noeaoouoa 

3.4  ft/asc  for  i?itd3 
reoca&aos 


%*  8  t% 


-  5S  * 

ft/oec 

V 

) 

.  juf  ,T._r.T,...r;. 

a-hji 

KfS 

Oil 

-.  4.$& 

i.v 

■  4.1 

0.? 

4.4? 

1.? 

ui 

1.3 

4.00 

2.0 

1» ,;) 

1.9 

3-53 

2.0 

1.4 

2.5' 

3.0S 

2.1 

1+  x 

1.1 

2.59 

2.5 

Iu4 

.3«r 

2.12 

2.5 

X«“4 

4.2 

1-73 

2.5 

1*5 

4.8 

1.25 

2.6 

\.  '”5 

Coeff ideate  (sete  pg. 

16  Of 

V  *  2.9  ft  /sec  for  heaire 

■  ras-jnBiieii- 

V  «  4.1  ft/sec  fox  pit.* 

rtesoaaacB' 


$»*  .follcsdag  v*iu*>»i  asssd$d  iTcs*  tfes  QGi%aWt&cm  <Sf  Is  sal  &*  urns 
uissfis 

it  a  1..8  St  (a  Xittlt  tw&jw?  &sa  lisOf  t§»  Usages,  a?  tlse  •v“»»3'8.\ 


1 


cSa  to  -ifee  yJ^SGfe*'***  -is 


ft,  ss  .25-'  ft  (sMtea  csf  ”0*$?  :j,!  sx-v^um^} 

fij  •"  .25  ft 

b  ss  ,.7§  ft  (dtfptfc  ax  g£3fcm*®Hwe) 

«3eiw5?> 

/  «* 

#ok-  different  ?»aeed8  V,  fte  \*iu«8  v  _,y£.  aj  and  *  y«-*~  t  «b£<&  tare 


.MHB 


el;*  madid  ea  iojpst  to  tSt^ocd^uttf  srwgsf-asj  aeye  eslcolwtad. 


V 

1.0 

2,0 

3.0 

4.0- 

5«2 

;  0.131 

0.263 

o.a* 

o,  526 

o.'635 

* 

I.3H 

.941 

•  >■? 

.20 

-0.26? 

(4  ”  4  ft) 

,73" 

,  0*999 

Mk 

.  .,6t9 

.444 

,.,222 

s 


Fig.  2  show*  the  result*  for  b  lad  3.  Values  at  b  and  B  were 
read  from  the  curves  for  the  various  speeds  that  were  selected 
in  the  Iterative  computer  program  to  find  the  motions.  However,  for 
the  vessel  in  8»ft  waves,  the  curves  showed  such  good  resemblance  to 
second  and  third  order  parabolas  that  curves  of  the  type  b»  qV  rani 
B  *  p'v3  +  <y'v»-  +  r'V  f  a*  ware  fitted  through  «Uffareat  points  of 
the  functions  b(v)  susd  B(v)«  The  result  uaa:  b  =  (0.0x55  V^- 
0.063  V  +  0. .152 )p  ,  B  -  (-0.00014?  V3  4-  O.lS  V2  -  0.034  V  + 
o.oa47),o. 

e  end  B  -  Coupled  Damping  Coefficients  (pg.  18) 

1=4  ft 

e  =  -0.3  V  +  0.01  V2 

E  -  -0.1  V 

1  -  6  ft 

e  =  0.3  V  -  0.01  V2  and  E  -  -0.1  V  f or  V  <  2  fii/sec 
e  -  -0.3  V+O.Gi  V2  and  E  s,  o.l  ?  for  V>2  ft/aec 
Dhoso  values  aro  very  disputable  since  no  theory  exists  nor  have 
experit-ar.tti  been  carried  out  to  evaluato  e  uad  E  at  sUb merged 
bodies  dose  to  the  surface  .and  oscillating  with  frequencies  near 
to  zero.  However,  the  e  and  E  values  of  the  basic  bull  arc  very 
small  in  comparison  with  the  e_  mid  caused  by  the  presence  of 
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A  «.  35*4  lbs 
I  -  0.0435  ft1* 

C  a  A  x  s  64  x  0.0435  «•  2.73  ft  lbs/rsd 

Fig.  4,  which  was  taken  from  Ref  .3,  shove  how  the  trim  changes  with 
forward  velocity,  the  curves  far  the  hydrodynamic  and  hydrostatic 
moments  at  a  particular  speed  are  drawn  in  the  adjacent  sketch  with 
the  assumption  that  the  hydrodynamic  moment  Is  a  second-degree 


V  >  °>>  /»„ 

Cte  the  basis  at  the  experimental  data 
of  Ref*  3  shown  an  the  adjacent  sketch, 
an  expression  relating  C3  to  speed  can 
be  found. 


If  we  put  C3(v)  <•  pv2  +  qV  +  rff,q,  and  r  can  ba  defined  by  the  aet: 


C3^vl)  -  4  *V1  +  r 

C^(V2)  •  p«f  4  qv2  ♦  r 

°3<V  *  P^  ‘  *v3  +  r 

the  result  la  ci(V)  -  -0.0917  V2  4  0.692  V  4  0.492.  The  expression  for 

o  «  _  « 

then  becomes 1  (0.0917  V2  -  0.692  V  -  0.492)^g*  *  c3  * 
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-e 

We  write  —  Instead  atJW since  the  hydrodynamic  moment  has  to 
have  the  opposite  sign  of  6  and  to  avoid  the  obvious  difficulty 
of  a  negative  6  under  the  square  root  sign.  If  the  values  of  Cj 
found  in  the  previous  paragraph  are  utilized  to  compute  the 
actions  of  the  semisubmarine  in  waves,  they  are  much  larger  than 
the  experimental  motions  of  Ref.  3.  Therefore,  apparently,  the 
values  Cj  determined  from  still-water  experiments  are  not  valid 
in  waves.  The  value  chosen  for  was,  therefore,  reduced  to 
-0.5  V  +  0.05  V2,  which  resulted  in  computed  motions  that  correspond 
fairly  well  with  the  experimental  results. 

5.  d  and  D  -  Coupled  Virtual  Mass  and  Moment  of  Inertia  Coefficients 

L 

?  0 

rrj  rZx  dx  was  calculated  according  to  Simpson ' s  Rule  and  was 
-L 

T 

found  to  be  equal  to  -0.0251  f ,  a  very  small  value, 
d  «  D  »  0.0251/)  . 

6.  g  and  0  -  Coupled  Spring  Constants,  pg.  21 
V  a  0.48  ft3 

m1  *  0.015 
D'  =  .001* 


0.152  V* 

Therefore,  g  ■  +  0.152  V2  lbs/rad 

i-8* 

X  a  — —  ft 
•  12 

Ay  *  0.24  ft2 

Mj  -  -1.18^  t  therefore,  0  =  +  l.lflfft  Ib/ft 


*  Based  on  deeply  submerged  submarine  data  ignoring  free  surface  and 
frequency  dependency. 
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7.  fluneriool  reeulte  for  pitch  cod  bunt  soplitudae 


X  *  4  ft 

V(f t/sec ) 

*1 

z2 

8^(  radians) 

®2 

0.1 

o.oo4 

0.004 

0.027 

0.024 

0.7 

0.009 

0.009 

0.033 

0.035 

13 

0.016 

0.016 

0.043 

0.037 

1-9 

0.098 

0.040 

0.056 

0.045 

2.5 

0.19 

0.21 

0.105 

0.10 

3-1 

0.125 

0..16 

0.15 

0.12 

3.7 

0.16 

0.32 

G.43 

0.35 

4.2 

0.25 

0.4 

1.2 

1.3* 

4.8 

0.21 

0.3 

1.0 

0.9* 

*  Unstable 


zj_  is  the  heave  amplitude  according  to  the  solution  of  the 
linearised  equation* .  sg  is  the  heave  from  the  timo-Motory 
solution  of  the  actual  differential  equations.  &  lo  the 
pitching  amplitude. 

1  -  8  ft 


V 

*L 

*2 

&i 

02 

0.1 

0.05 

0.1 

0.084 

0.09 

0.7 

0.07 

0.13 

0.10 

0.10 

1.3 

0.075 

0.125 

0.11 

0.10 

1.9 

0.112 

0.17 

0.145 

0.13 

2.5 

0.35 

0.35 

0.8* 

e.26 

0.2 

3.1 

0.29 

0.27 

0.3 

3.7 

0.17 

0.5 

0.39 

0.5 

4.2 

0.22 

0.38 

0.7 

1.0* 

4.8 

*  Utaetdbla 

0.56 

0.6* 

2.0 

1.1* 
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AfflfflDIX  V.  EFFECT  CP  reiCTZOB  CP  THE  M0IX0«3  CP  A  SOC80BWBIXS 
A a  noted  on  pg.  6  of  the  text 

Cfrlction  baa  little  effect  on  the  motion  at  a  ship  in  a  seway. 
Ihere  if  a  frictional  resistance  proportional  with  the  square  of  the 
velocity  In  the  direction  of  the  notion.  This  increase  in  frictional 
resistance  is,  however,  very  snail  with  respect  to  change  Id.  flow  a 
periodically  saving  Alp  creates  around  it. 

For  a  eanl  submarine  it  could  he  a  different  natter  since  the 
sail  can  have  a  relatively  large  area  and  its  drag  applies  at  a 
considerable  distance  fron  the  center  at  buoyancy.  If  the  dreg  is 
represented  by:  0  *  1/2  p  V2  Cj  S  where  S  la  thn  underwater  area  of 
the  strut,  Cp  has  an  average  value  of  O.Ol  for  reasonable  roughness 
of  the  strut.  This  gives  for  V=  k  ft/aecj  0  =  1/2  x  1.99  x  16 
x  0.01  x  2  x  0.5.  D=  0.16  lbs.  This  la  for  n  Froude  Humber 
O.35.  The  pitching  oaaast  la  then  0  x  0.5  *  0.00  ft  lbs  which  is 
negligible  in  comparison  with  the  moment  due  to  th»  flow  around  the 
body.  The  hydrodynamic  derivatives  that  are  influenced  by  this  drag 
are  5^,  1^,  X,,  M,,  and  They  can  be  evaluated  as  followst 

*U=  -P  «(Cd8  +  C*  Sb) 

l^=  -p  U  Cjj  3  h, 

2u-  p  U  CLb  St 

Xg  =  1/2  p  Cp  1#  V2 

\  =  l/Sp^^/.pgx^ 
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ttea  Cj^  s  drag  coefficient  of  the  body  of  revolution 

Cj^j  -  lift  coefficient  cf  the  body 

-  equivalent  hydrofoil  area  of  the  body 

hg  =  height  of  point  cf  application  of  drag  on  atrut  above 
the  center  of  buoyancy 

lg  -  lez^th  of  atrut. 

An  evaluation  of  these  coefficients  item  that  they  care  email  in 

due 

comparison  with  tha  other  derivative*!.  In  the  pert^to  the  drag  cf 
the  strut  will  be  much  smaller  then  p  g  x8  Ay.  2a  first-order 
approximation  Xy,  Z^,  and  Xs  nay  be  neglected. 
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-  Agggjgn  vx.  on  cowhol  sraxtots 

Hie  principal  requirementa  for  stabilizing  fine  are  symmetry  of 

the  croea  section,  high  nsuclnxun  lift,  lew  dreg,  resistance  to 

cavitation  and  great  strength.  It  follows  frota  the  latter  consideration 

that  tbs  plan-form  must  be  fairly  stubby,  with  an  aspect  ratio  usually  not 

larger  than  two.  This  low-aspect  ratio  has  also  the  advantage  of 

making  the  angle  of  stall  larger  although  no  increased  maximum  lift  is 
thereby  achieved. 

Lifting  lira  theory,  applied  to  a  syanatrical  hydrofoil  of  aspect 
ratio  /?>  gives  for  the  lift  coefficient: 

°L  »  2n<*:  *  ■&>• 


»cL  1 

Thus,  — —  3.  .  'jfrjia  theory  oeouEea,  however,  that  the 

a  *  -±-  +  -i- 

2t»  &T 


lift  distribution  is  elliptical  and  that  the  plan  form  is  also  elliptical. 


The  theory  of  lifting  surfaces  brings  in  a  aarroctioa  for  rectangular  wings 


“4  for  U»  uglc.  of  Mttek. 


T 


sin  or.  +  2  oin  <c  where 


i  Be 

2ir  tt-/?. 


E  »  1  + 


"aF  ' 


Zf  the  aspect  ratio  A  is  taken  to  be  equal  to  2,  vs  get 

acT 


B  m  1.4o.  C.  a  2.6  sin  o  +•  2  sin*  o  cos  oc,.  *  2.6  cos*1-  +  4  (sin  cc 

d  a 

cos  «)-  6  sin^  oc  a  +  2.6  +  4  cs.  -  6  if  sin  oc  is  approximated  by  ac  and 

cos  «  by  1.  A  free  surface  correction  is  not  necessary  since  the  foils 

are  sufficiently  deeply  submerged. 

jSk.,.1. 6  sin  a.  +•  4  cos  a.-22  sin2  «  cos  cx 
»  a?  8 

»  -  2.6  oc  *  4  -  22  ct. 


-  6o 


If  we  expand  Cg,  into  e  Taylor's  aeries  up  to  the  second  order,  we  get: 

3c.  d2c.  9  . 

Ct  «  - -  «  +  1/2  - a2  *  2.6  a  +  6  «  -  1.3  a?  -  17  *  . 

da  da 

This  lift  coefficient  takes  on  the  following  values  for  different  a' at 


a  (radians) 

a0 

°L 

cv« 

0.1 

5°Uo* 

0.32 

3.2 

0.2 

n°30' 

0.72 

3.6 

0.3 

17°20' 

1.16 

3.85 

If  we  assume  the  foils  to  operate  alauaoifiaHy,  the  lift  -on  a  foil  cca  he 
expressed  as  Lq  ain  o>et  since  the  lift  is  proportional  to  the  harmonically 

changing  angle  uf  attack. 


ac. 

To  find  an  average  — — — ,  which  will  causa  the  foil  to  operate  in  euch 
3  a 

s  way  that  the  work  done  during  a  quarter  of  a  cycle  will  be  the  sane  ao 


with  a  changing 


one  can  reason  as  follows.  Assume  the  ship  to  move 


3  a 

harmonically,  the  pitch  angle  can  he  written:  6  =  0Q  cos  <e  t,  0  *  -  60  u>  sin  <»t. 
If  the  anti 'pitching  control  Is  of  the  form  a  *  k^e,  we  get  : 


«  *  V  0a> .  sin  O  t  a  sin  a>  t- 
The  work  done  by  the  foil  during  a  quarter  cycle  is: 

, .  rTA  . 

W-J  L9dt 

j0  /T^  x 

=  J  pAV  CL  oCa  si*  cut  dt 

°  st/n 

-=.  Const  \  CL  Sit  U)t  cOz 


With  Cr  =(  I.  .^a  where  is  taken  to  be  constant,  w  geti 

^  'J  a'»  Vf 

W'  =  C.,>t  5i*  *"  £ot  cit 

~c°"‘ 

(ich  \ 

If  we  put  V  -  W,  w  can  find  the  unknown/— —  I 

l  3  ct  /*• 


n \ 


-  - 


J  (2X0I+  Ux~  UU3-  17^  ^  <<t  dt 


u  JE 


4w 


/  acL\ 

When  we  substitute  «  ~  or.  ain  o>  t  in  this  expression,  (  -  ■  -)  bscotaea  equal 

3  «c  M 

to  3.73.  For  a  real  fluid,  however,  the  lift  coefficient  will  be  lower. 

The  lift  needs  a  certain  time  to  build  up  completely  so  that,  since  the 

foil  oscillates  constantly,  there  will  be  a  loss  in  lift.  Cta  the  other 

hand,  separating  of  the  boundary  layer  also  needs  tied  so  that  the  angle  of 

stall  la  higher  with  oscillating  fins. 

Considering  these  end  other  unsteady  effects,  it  la  reasonable  to  assume 
3Ct 


tbs  equivalent  a  3.U,  Ibis  value  wee  chosen  to  compute  the  idealized 

3  ot 

control  parameters  *0,  l±,  t,  ^  (000  Chapter  IV)  and  in  the  first  trials  of 
the  more  realistic  control  system.  The  urea  of  the  fins  was  taken  as  2.5  per 
cent  of  the  longitudinal  cross- oectional  area.  Five  per  cent  le  a  maximum 
value  far  rigid  anti -pitching  fins  in  surface  ships  (Fof.  18).  It  was  found 
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from  tbs  computational  results  that  a  sranitubesrine  with  foils  with  the 

previously  mentioned  characteristic*  sod  dimensions  can  ho  stabilised  la 

U-ft  wavs*  (X  •  L)  hut  not  in  wavs*  of  double  that  length,  (see  pg.  37). 

As  a  first  naans  of  increasing  ths  lift,  the  area  of  the  fins  was 

increased  by  30  par  cent.  Ohio  is  allowable  in  this  ca^s^noe^lMf*8 

pitching  moments  are  much  smeller  and  act  with  lower  frequency,  tbun  in  ahead 

waves  where  the  frequencies  of  encounter  ere  higher  end  the  moments  ere 
In  ahead  seas 

higher,  fit  will  he  necessary  either  to  retract  the  folia  or  activate  than 

such  that  the  lift  remains  eero.  It  must  he  noted  that  no  stabilisation 

is  necessary  on  the  semiauhoarlno  in  ahead  -oaves, (Ref.  1  and  3). 

A  second  means  to  improve  the  performance  of  the  foil  consists  of  the 

uee  of  flaps  which  increase  the  lift  coefficient  greatly.  Data  from  Ref.  (17) 

show  that  it  is  possible  to  obtain  =  1.5  at  18°  tilt  where  both  the  main 

part  of  the  fin  and  the  flap  are  activated. 

Cavitation  can  become  an  Important  factor  at  ship  speeds  above  25  knots. 

The  lift  coefficient  needs,  therefore,  to  be  decreased  sufficiently  (Ref.  19). 

9Cl 

A  CL  of  1.25;  ■■■■*■  **t  seems  a  reasonable  assumption.  This  value  lies 

a  little  above  the  CL  found  in  teats  conducted  at  Davidson  laboratory  on  a 
model  with  oscillating  anti-pi tching  fins  (Ref.  20). 
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AgPSHDDC  VII.  XBBAt  CCMIBCE, 

(see  pg.  30~-  33) _ 

tii":  wlnee  s*  iti-  «ad  m-  gvvea  in  A&5*  axprefi^-ian 


of  Chapter  IV ,  pg.  31>  becomes*  for  X  ■  4  ft: 

O  (-0.642  +  0.0008  V)d  +  (“0.9  +  0.07  v  *  0‘5®  h  v) 


46.1  A.2  V4 

<SS»  longitudinal  cross- eectiaonl  wwa  is  1.2  ®1  so  t»t  if  vs  A 
eqjjfel  to  5  jar  cast  of  this  area  (Baf.  2B),  A  teacssss  O.Ov.  s<|  ft.  Si»U8f«> 

O  (-0.042  +  O.CCOS  V}’-  1-  (-0.9  +  0.039  vr 


A  £5C  —  . 

o 


0.166  v 


p 

For  V  *  4  ft/aso,  6q  «  . 01335 » 

4  *  0.115  *  5°40'. 
o 

2 

V  n  8  ft/mu,  60  -  .259, 

«s.  -•  C:.  ‘>06  =  29 J 

V/ 

In  8-ft  waves  (1  »  2L),  ws  get : 


2  (0.0§4  -  0.02  V  -  1.22  A  vf  +  (I.538  +  .15  V)2 

6o  *  - - - — - - — - — 

46.1  A2  v4 

If  we  take  the  same  values  for  A,  the  following  results  are  obtained: 
V  »  4  ft/sec,  «  0.109 

io  »  .33  a  19° 
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V  *  2  ft/eec,  82  =  1.275 

6q  *  1.13  radians  »  55° 

It  lo  seen  that  complete  pitch  stabilization  Is  theoretically 
possible  in  4- ft  waves  al though  the  fin  angles  have  to  become  rather 
large  with  low  forward  speeds.  In  waves  of  double  the  Ship  l»J«th, 

complete  control  become  harder  at  the  high  speeds  and  is 

impossible  at  lav  speeds  (Fr  »  .16). 

Far  complete  heaving  control,  expression  (3),  pg.  33,  for  the 
motion  in  4-ft  waves  gives: 

2  (0.15  +  o.o4  v)*  +  (o.48  +  0.016  vf 

*i - - - 

0.0388  vk 

where  the  foils  are  the  same  as  in  the  previous  case. 

For  V  *  k  ft/sec,  8^  a  0.04l6 

^  a  0.204  *  12°. 

At  the  speed  of  2  ft/sec,  8^  *  0.504  <*  29°. 

In  8-ft  waves  the  expression  becomes: 

a  (-1.147  +  0.23  V)2  +  (0.87  -  0.012  V)2 

8^  a  'miM  . . mi  . .  . . . . . . . . .  m.—i  •'■■■ 

0.0386  V4 

Fear  V  *  4  ft/sec,  82  «  0.0726 

*  0.27  =  15. 5°- 


reap  V  *  2  ft/sec,  6*  «  1.92 

*  1.385  *  79-5°. 

Again,  It  1«  impossible  to  control  the  notions  completely  at 
speeds  in  long  w%t. 

Zt  should  be  noted  that  the  foils  used  in  these 

calculations  ere  very  large  but  that,  on  the  other  hand,  the 
could  be  increased  by  the  use  of  flaps  (App.  VI). 
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AFPKHDIX  Vin.  COWER  OL,  gfSCSIEJ* 

In  Chapter  IV  (pg  34)  tbe  following  expressions  m  attained: 

e,  »  1/4  p  tt(U  +  P) 

e  l/  4  p  t?C  ^  ^ 

=  l/4pn(Hscg2  +  fccj^) 

In  a  rirst-trial  a  foil  area  of  2.5  per  cast  of  the  longitudinal  arose- 
sectional  area  was  chosen,  sg  eg  «  Bp  cp  »  0.03  ay  ft,  where  eg  la  span 
and  eg  la  chord  of  tbe  anti-beaaring  foil  end  analogously  for  the  antl- 
pitchimg  foil.  With  an  aspect  ratio  of  2,  we  get: 

0j  =  Sp^'  0.245  ft 

CH  =  CP  ~  °-122k  ft 

so  that  a.  =  0.012 

f 

Af  *  0.024 

d.„  =  Dp  «  o.OlB 

If  we  use  similar  reasoning  as  fear  hull  on  pg.  15,  l.e.,  that  with 

decreasing  frequency  at  encounter,  the  virtual  mass  and  mo— nt  of 
Inertia  Increase,  we  am  writer 

a  e  0.012  +  C.003  V 

f 

Af  a  0.064  ♦  0.006  V 

With  this  assumption  the  Virtual  case  end  moment  are  about  doubled  at 
the  speed  of  aero  frequency  of  encounter  compared  to  zero  speed. 

bf  «  2  x  1/apA  V/— JL  *■  ^ 


«  0.204  V  If  — i-  taken  to  be  3-1*  and  CB  »  0.02 
&  C 


*  A  summary  of  the  different  control  systems  mentioned  in  this  Appendix 
can  be  found  on  pg.  75»- 
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•f  *  1/2 p A  L2  v(— £-  f  cX  ,406  V 

«.S*  «l/l»p*VL  f— L  ♦  Ctj)  tt  .204  V 
r  *  '  &  «  7 


z4  *  l/2pA  V2  +  Cj^-e  0.103  Vs 


Mg  -  0.206  V2 

A  tins  lag  of  0.1  adt^wds^  taken  into  account  for  the  four  control, 
elements.  Hie  method  of  solving  the  linearieed  differential  equation!) 
was  a  trial  and  error  method  in  which  four  values  of  the  different 
control  parameters  were  taken.  Hie  coefficients  cf  the  six  per  asm  tent 
i,  *,  1,  %,  6,  6,  were  computed  and  Hurvitx's  Stability  Criterion 
applied; 

Hw  two  motion  equations  without  exciting  forces  or  momenta  can 
be  written  mm 

[/tn.+o.)®Z+  J&D  +  cjx  +(i||)l+  =0* 

[D©t^EO  +  Gjx4-[(j-tA)(D,'+6(D<f<!j0  *  o  * 

u»  characteristic  determinant  of  uus  set  is  of  the  farm: 


Vi  +  Vi  f  Vi  +  +  3  0 


11 


His  neoessery  and  sufficient  conditions  that  the  motions  will  he  stable 
is  that  all  cr^'s  be  negative.  This  will  he  ao  if 
ai  »  0  (1  =  0,  ...,  4) 


®1  *0 

*3  "a 


>  o 


and. 


*1  *0  0 

®3  *2  *1 

0  *4  *3 

d2 


*  ©  stands  for  jT  and©  for 


>  0 
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Dm  values  at  the  control  parameters  that  caused  the  coefficient* 

a,  b,  c,  A,  B,  C,  ...  to  be  such  a*  to  fulfill  these  condition*, 

war*  used  to  coogute  the  pit  Thing  end  heaving:  motion*.  Due  optima 
.  already  been 

value*  have^^nted  out  sed  showed  the  follavitc  result*  in  4-ft  waive  i 


(Linear  Solution)  Table  1 


vfcy**c) 

h 

*2 

*3 

K 

afre) 

e(rad. ) 

2-5 

-0.7 

•0.8 

0.95 

-0.58 

0.1 

0.035 

3*1 

-0.7 

-1.92 

0 

-0.37 

0.08 

0.09 

3.? 

-1.4 

0.3 

0.24 

-0.2 

0.076 

0.11 

4.2 

-2.6 

-1.7 

0.7 

-3.9 

0.072 

0.17 

4.8 

f  -2.9 

-1.4 

0.5 

-0.9 

0.013 

0.13 

Aden's  aumeric&l  method  of  solution,  which  take*  the  nonlinear  effects 
into  account,  gam  the  following  results: 


Table  II 


V 

*1 

*2 

*3 

k4 

t 

fcr 

2.5 

-1.3 

•0.4 

0.1 

-2.0 

0.10 

0.08 

3.1 

-1.4 

1.8 

1.0 

-1.2 

0.12 

0.10 

3.7 

-1.8 

-0.3 

-0.1 

-1.0 

0.12 

0.16 

4.2 

-3-5 

-1.0 

0.9 

-6.0 

0.15 

0.15 

4.8 

-5.1 

0 

0 

—2*0 

0.33 

0.12 

With  a  time  lag  twice  ea  large  and  a  foil  area  3.5  par  cent  of  the  cross- 

sectional  area,  we  get  with  the  same  method: 

Table  III 

V 

*1 

*2  *3 

*4 

* 

6 

2.5 

•1.2 

0.4  0.2 

-2.0 

0.02 

0.033 

3.1 

-1.3 

-1.4  1.0 

-1.2 

0.004 

0.042 

3.7 

-1.7 

•0*3  0*1 

-1.0 

0.016 

0.06 

4.2 

unstable 

4.8 

unstable 
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At  low  ipMdi  the  larger  foils  are  more  effective  the 
foils,  but  at  higher  speeds  the  tine  lag  has  such  a  large  influence 
that  the  virtual  mass  and  moment  of  Inertia  become  negative. 

In  6-ft  waves  the  linear  nethod  showed  the  foilowii«  results: 


Table  iv 


V 

*2 

*3 

z 

6 

1.9 

-1.0 

«*0s  4 

0.5 

-0.5 

0.08 

0.15 

2.5 

-l.l 

-0.4 

0.5 

-0.4 

0.05 

0.16 

3-1 

-l.o 

0.29 

-0.13 

-0.1 

0.13 

0.18 

3-7 

unstable 

4.2 

-1.64 

-0.4 

0.4 

*1»  4 

0.10 

0.20 

4.0 

-2.2 

-0.1 

0.1 

-0.7 

0.11 

0.18 

Here  a 

time  lag  of  0.1  sec  was  assumed  and  the  area  of  the  foils  2.5 

per  cent.  The 

(Adam's  nume 
second  j^cethod 

vaaTi^tterded  into  confuting  the  angles  cf 

attack  on  the  foils  in  order 

radians 

to  check  if  the  angles  do  not  exceed  a 

may i man  of  0.35  to  0.4jjAich  is  allowable  for  dimple  flna  of  snail  aspect 

ratio. 

The  following  results  were  obtained: 

Table  V 

V 

kl 

*2  *3 

fc4 

z  e 

«p(rad)  8^ red) 

1.9 

•1.0 

-0*3  -0.3 

-1.0 

0.07  0.17 

0.4 

0.2 

2-5 

-1.0 

-0.3  -0.3 

-0.9 

0.15  0.20 

0.4 

0.3 

3.1 

•1.2 

-0.3  -0.3 

"1.0 

0.28  0.25 

0.4 

0.4 

3.7 

-1.2 

-0.4  -0.4 

-1.2 

0.26  0.35 

0.56 

0.2 

4.2 

•1.6 

-0.4  0.8 

-1.5 

0.06  0.25 

0.5 

0.35 

4.8 

-2.2 

-0.1  0.1 

-2.2 

0.09  0.22 

0.45 

0.2 
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6^  :i»  the  angle  of  attack  ;jn  the  anti-pitching  fin,  not  taking  into 
account  the  wave  particle  velocity.  Analogously,  for  ftg.  Ibe  matters 
In  these  columns  are  the  nsKiEua  angles  of  attack  the  foils  take, 
according  to  a  linear  O^Jn,  In  reality,  there  will  be  overloading 
when  or  &H  >  0.4.  Analogously,  with  ths  aotioan  in  X  «  4-ft  wave*, 
the  festive  and  pitch  beear©  unstable  at  speeds  higher  than  3  ft/esc  when 
the  lag  Is  doubled. 

The  vessel  with  two  anti -pitching  flna  of  1.75  per  cant  area  each 
mud  cm  pair  of  anti -heaving  flna  in  the  middle  of  total  area  3.5  par 
cent  has  the  following  control  characters 

■2f  *  0.019  +  0.004  V 

Af  *  0.036  4  0.008  V 


a  0 

Cf  -  0  With  tale-  t?rn,  whicl  war  ;;  ositiv*' 
xt.  Xhx  previous  ease,  v?  lore  a 
certain  amount  of  our  ritcai.3'-’ 

?»•»  *ye  ■“  -1-1  -V  -}•  fT  -v^r>, 

Zg  *  0.142  V2 

«  .20  V2 


It  vos  ato'5'loa®  Xvcm  the  ypwiaaa  reyaita  that  tite®  %  s»  gain  in 
having.  ths  haava  also  usstrol  tfc®  angle  of  tste  aatl^ltcSsiag;  foil,  dost 
the  pitch  emtrol  the*  eogl©  <£  tha  SBti-tozawiag  foil.  SsNfoni  the 


centre!  was  stagUfied  to: 


5  =  1^6 
1?  i 


!s,t8 

cj,  A 


8-  =  lv.,K  -  &  t.  Si 

a  *?■  <t  X 

where  t^.  vee  tshen  to  he  0.1  ace.  Sse  ifoUoeSca  results  vs»®  dhtaHseds 

Table  VI 


V 

kl 

St. 

tj. 

*/ 

-1 

V 

f* 

r 

1.9 

-1.1 

-1.1 

0.02 

0.10 

0.3? 

0.18 

2.5 

-1.1 

-1.1 

0.12 

0  * '  ' 

0. 4? 

0.43 

3.1 

-1.8 

-1.8 

0.2 

0.15 

0.42 

u.tkj 

3.7 

-1.3 

-1.3 

O.lQ 

0.17 

0.51 

0.46 

4.2 

-1.4 

•=1.4 

0.2 

0,21 

0.46 

o.  ko 

4.9 

-1.5 

-1.5 

0.3 

0.26 

0.43 

0.49 

«»aj  r>vaiif«.  4*i(<.r  /WKira 


The  values  for  <o  were  toners,  aa  usual,  ae  the  aastoa  'value  of  pitch- 
angle  aut-to-out  divided  by  two.  The  average  pitch  eagle  erouad  which 
the  vessel  was  oscillating  had.,  however,  the  tendency  to  bocoaa  larger 


end  larger. 

dynamic  pitching 

To  eliminate  the  loss  instability,  flBp~fina  were  introduced, 

while  the  total  fin  area  was  increased  so  that  the  gain  in  control  io 

dCi 

doubly  increased.  The  — •  van  aesuoud  to  be  4.8  (*i«e  App.  VI),  the 


9  m 

area  of  the  pair  of  anti -heaving  foils  3-5  per  cent,  the  area  of  each 
anti -pitching  foil  2-5  per  cent.  Furthermore,  a  step*  control  was 
assumed  (aeo  Chapter  IV).  The  result  is  given  on  the  following  page. 
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Table  VII 


V 

*1 

*4 

z 

6 

1 

fF 

1.9 

-0.3 

-0.3 

o.o6 

0.02 

0.08 

0.03 

2-5 

-0.3 

-0.3 

0.05 

0.006 

0.07 

0.07 

3.1 

-0.3 

-0.4 

0.001 

0.005 

0.05 

0.03 

3*7 

-0.3 

-0.4 

0.001 

0.006 

0.05 

0.05 

4.2 

-0.3 

-0.4 

0.002 

0.008 

0.04 

0.05 

4.8 

-0.35 

•0.4 

0.002 

0.003 

0,03 

,0.04 

This  control  appears  to  be  too  "stiff."  The  accelerations  amount  to 
twice  and  even  three  time a  the  limit  of  passenger  comfort.  The  6^  and 
?B  indicate  the  actual  angle  of  attach  between  the  foil  mats  stw  the. 
water  particle  Whexas  ami  4^  of  the  pravi«jus  re  suite  did  not  tabs  the 
direction  of  the  wave  particle  Into  account. 

Again,  a  trial  van  cede  with  one  entlrplttiblng  foil  (at  the  stem) 
and  a  pair  of  anti-heaving  folio  midship.  The  eras  of  tfea  latter  wen 
not  changed  (3*5  per  cent),  while  the  area  off  tins  foraer  vae  inoreatae 
to  3-5  per  oent. 


Table  VIII 


V 

*1 

*4 

S 

6 

*; 

1.9 

-0.3 

•0.4 

0.09 

0.08 

0.36 

0.66 

2.5 

-0.3 

-0.4 

o.i4 

0.07 

c.35 

0.60 

3.1 

-0.3 

•0.4 

0.19 

0.065 

0.28 

0.60 

3.7 

-0.3 

•0.4 

0.25 

0.16 

0.33 

0.62 

4.2 

-0.2 

-0.4 

0.10 

0.11 

0.28 

0.48 

4.8 

-0.2 

-0.4 

0.004 

0.063 

0.25 

0.48 

The  veeeel  with  two  ecall  anti-pltching  fins  (1.6  per  cent  area  each, 

— —  -t  4  ),  directly  attached  to  the  body  of  revolution,  has  the 
d  cf 


following  charaaterietice: 


-  73 


Sf  =•  0.009  4  0.004  V 
kf  a  0.038  ♦  0.006  V 

df  a  s  0 

®f  -  Bf  »  0 

Cf  a  0 

bf  a  .29  V 
Bf  a  .68  V 
Zg  a  .120  Vs 
>%  a  .346  v2 


A  etepcuntrol  by  9  and  s  «u  applied,  end  the  following 
amplitudes  and  parameters  were  obtained i 


Table  IX 


V 

kl 

k4 

z 

e 

■a 

*E 

1.9 

-0.3 

-0.4 

0.065 

0.063 

0.27 

0.6 

2.5 

-0.3 

-0.4 

0.16 

0.1 

0.29 

0.53 

3.1 

-0.3 

-0.4 

0.25 

0.14 

0-3 

Oe  ^v2 

3.7 

-0.3 

-0.4 

0.13 

0.10 

0.29 

0.52 

4.2 

-0.3 

•0.4 

0.06 

0.09 

0.32 

0.47 

For  the  high  speeds  the  pitch  sod  heave  still  have  the  tendency  to 

oscillate  around  an  increasing  positive  or  negative  value.  Therefore, 

in  the  next  try,  a  supplementary  position  control  vas  introduced. 

$  z 

*P  ”  kl  75/  +  «H  3  k3z  +  k4  727  * 
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Table  X 


V 

*1 

k2 

> 

k. 

4 

j  *  • 

e*  -<S)  A  f 

|  F  H 

1.9 

1 

0.06 

0.096 

2.5 

I 

0.27 

0.096 

3-1 

-0.S5 

-1.2 

-1.2 

-0.25 

0.30 

o.o46  *e  0.3 

3.7 

0,17 

0.090  j 

4.2 

0.09 

0.060  ! 

4.3 

0.02 

0.020  ! 

! 

Use  6  sanies  occur  only  at  -she  asoaent  v>Sieia  the  foil  is  flipped  over 

H 

fraa  am  wsxtnum  to  t be  other.  Since,  In  pgectice,  tfc®  foil  needs  a 
certain  time  to  flip  over,  tie  actual  aastasa  aagSa  of  attack  will  be  loos. 
This  system  of  control  applied  m  the  vassal  moving  la  Wa  waved  AoweS 
pi tcli  amplitudes  Usat  vsre  all  bale-;?  0,05  ®i*3  heave  amplitudes  below  0.10. 
rMn  was  to  b®  expected  since  exciting  faros  and  moment  are  such  renal  lag 
than  in  the  8-ft  wave  cane.  Contrary  to  the  motion  of  surf  ecu  oJiipa  in 
abend  waves,  'share  wive®  of  length <?.mm  the  severest  aotiana,  the 
senieuteaaariaa  experiences  them  in  waves  of  about  twice  Ibe  reason 

for  this  Is,  besides  the  one  c&nticnod  before,  that  tbs  frequency  of 
encounter  Is  such  that  for  these  'saves  the  spaed  of  resonance  (Fr — -3)  lies 


*  75  “ 

in  tbs  range  of  acsemiL  *Mp  tasgar  ’aswea  csuj®  &i$s@r  sgnsafl®  o£’ 

resonance  atoll*  aitfe  stoartesr  acres  t&a  jfeeoas'a®;-  siysea  la  tmss. 

Fig.  5"  dboae  tisa  results  of  ttoo  nmt  control  eyatens 

as  cctt$fiore<l  with  the  a»3Sontrel3is4  mtlsma. 
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appskpu  ix.  ctounaa  frogbamb 

it  would  bm  iapoaaiole  to  give  is  justing  or  all  tan  program 
that  w*r«  written  to  canputo  the  amplitudes  of  the  motions  without 
and  with  control* 

A  set  at  four  exaqplea  ha*  b*en  chosen  to  give  an  idea  of  how 
the  heave  and  pitch  amplitude*  were  found = 

Ohe  first  example  computes  the  heave  and  pitch  for  9  different 
speeds  in  astern  waves  of  4  and  8  ft  length,  according  to  the 
fetniliar  method  at  harmonic  solutions  of  linearised  differential 
equations. 

Synhols 

A(l)  =  .a  +  a  3(1)  =  b 

A(2)  =  +  A  B(2)  it  B 

c(l)  a  C  D(i)  IS  6 

C(2)  =  C  D(2)  M  £ 

E(l)  *  a  0(1)  «  s 

B(2)  a  E  G{8)  *  G 

AF  e  coefficient  cf  cos  o>  t  in  expreeaioa  of 
total  exciting  fcrce®F 

BF  9  coefficient  of  sin  u^t  in  expression  of 
total  exciting  moment  M 

CM  »  coefficient  of  cos  u>  t  in  expression,  of  M 

DM  *  coefficient  of  sin  u^t  in  M 

THETA  (I)  *  angle  at  pitch  assumed  at  ship 
speed  V(Z)  (see  App.  IV) 

Z  s  heave  amplitude 

TBET  -  pitch  amplitude 

FIX  a  phess  leg  between  wave  end  heave 
aaxlaua> 


FUS  a  phase  lag  between  wave  and  pitch 
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In  the  second  asampln  tho  motions  sore  calculated  according  to 
Aden's  aurasrieal  Method.  On  subroutine  KESPOM  sets  the  initial 
conditions  (tbs  position  of  the  vessel  at  T  *  0).  It  ccaputas  F 
end  M  in  function  of  the  variable  T  sod  solves  the  two  motion  l 

equations  for  9  end  *.  ten,  different  library  functions  are  called 
(DPHVF  end  33S2VF)  which  cauputa  the  actual  %  md  &  nt  t iwt  time  T. 
Thereafter,  T  is  incra&sed  by  ’XL’;  and  the  sequence  etarta  all  over 
again. 

Tfco  third  program  is  eoasatiaily  the  sat*  m  the  first,  but  tho 
control  parameters  era  incorporated  in  the  coefficients.  Tfeo  two 
equations  are  then  subjected  to  Eouth-Ifeiralts's  Stability  Criterion. 

The  equations  are  then  solved  for  th®  values  of  the  control  parezustera 
that  fulfill  the  stability  criterion  in  the  sama  manner  as  the  first 
program.  The  control  psartvefisrs  are,  respectively,  PESO.,  PXSE2,  PKK3, 
fKXb.  To  each  of  ther**  values,  bJC's  ere  added  la  such  a  vay  that  the 
equations  ere  solved  for  ten  thousand  different  combinations  of  the 
control  parameters.  In  other  wards. 


First  Combination 

pna 

ms 

TO 

ms* 

Second  Ccssbtnatiaa 

Pisa  +  aa. 

mz 

MX3 

PICK-1* 

Third  Combination 

ma  *■  2iaa. 

psx2 

STO 

PK& 

Tenth  Combination 

FKKL  +  9DJQ. 

F WS 

TO 

FEE4 

Eleventh  Combination 
etc. 

PKKL 

PEE  +  £SC2 

SKK3 

m3* 

The  values  at  <z  and  9  that  are  below  the  specified  values  R(I) 
and  V(I),  respectively,  ere  pointed  oat,  togetlsar  with  their  pteee  1 ew® 


Xt  gives  a  tine  history  at  pitch,  heave,  sod  Xhts  respective  eagles  a t 
attack  on  anti -heaving  and  anti -pitching  foils  (D2LB  and  BELP). 
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♦HEAVE  AND  PITCH  WITHOUT  CONTROLS  LAMDA=4AND8FT » L I  NEAR IZED  EQUATIONS 

*  XEQ 

*  LI  ST  8 

C  HEAVING  AND  PITCHING  RESPONSES 

DIMENSION  V ( 30 ) *  A 1 1  30 > . A2 ( 30 ) .B 1 ( 30 ) , 82 ( 30 ) » A ( 2 ) *B ( 2 ) *C ( 3 ) »D ( 2 ) * 
2E ( 2 )  *G ( 2 )  , THET  A ( 30 ) 

PRINT  5 

5  FORMAT  (8H  LAMDA=4 ) 

READ  1 .  ( V<  I  )  ,B1  (  I  )  »B2U)  »A1  (  I)  »A2<  I  >  *1  =  1.18) 

READ  1.  (  THET  A  (  n.I=l,18) 

PRINT  1. (THETA! I ) ,1=1,18) 

PRINT  1 »< V(  I  )  ,1  =  1.18) 

1  FORMAT  (5F14.8) 

1  =  1 

C(1 1=7.74*1.99 
D(l) =0.0251*1.99 
D< 2 )=D( 1 ) 

G( 2 ) =2.43 
PRINT  2 

2  FORMAT  (75H  V(I)=  Z  THETA  F 

2LZ  FLTH  ) 

53  C<  3  > =0 • 09 17*V ( I ) **2-0 . 692*V ( I >-0.493 

C(  2  >  =  <2.78*THETA< I ) +C <  3 ) *SQRTF ( THET A (  I ) > ) /THETA ( I ) 

G ( 1 )  =0. 076*1. 99*V< I  1**2 
20  A ( 1 ) =  A 1  (  I  ) 

A< ? ) =A2 (I) 

103  I  F  (  1-9)70.70,71 

70  AF=0. 072-0. 00816*V( I ) 

BF=0.13+0 . 04* V ( I ) 

CM=-0.9+0.07*V ( I ) 

DM  =  -0. 042  +  0. 0008*V(  I  ) 

OMFGA=1.57*(V( I  1-4.54) 

B ( 1 ) =B 1 < I ) 

B ( 2 )=B2< I ) 

E(1 ) =-0 • 3*V (  I  ) +0. 0 1*V (  I  )  **  2 
E  <  2 )=-0.1*V(  I  ) 

GO  TO  26 

71  PRINT  102 

102  FORMAT  (8H  L AMDA=8 ) 

AF=0. 133-0. 0104*V(  I  ) 

BF=-1 • 1 +0 • 2*V ( I ) 

CM=-1. 538-0. 15*V(  I  ) 

DM=0. 094-0. 0206*V<  I  ) 

OMEGA* 1.5 7* (V( I  1-6.4) /2. 

B(1 )  =  ( 0 . 0 1 5  5*V ( I ) **2-0.08  3*V( I >+0.1521*1.99 

B ( 2 )=<-0.000147*V( I >**3+0.0 1 603*V ( I ) **2-0 . 0372*V ( I ) +0.0247 ) *1 .99 
E(1 >  =0 . 3* V ( I ) -0 . 0 1* V  < I  1**2 
E<  2  > =-0. 1*V (  I  ) 

*  SOLUTION  OF  SET  OF  DIFFERENTIAL  EQUATIONS 

26  P 1 = - A ( 1 ) *OMEGA**2+C ( 1 ) 

XF  =  ATANF<  DM/CM) 

XX=CM**2+DM**2 

SMl=SQRTF(XX)*COSF(XF) 

SM2=-SQRTF(XX)*SINF(XF> 

YF=ATANF(BF/AF) 

YY  =  AF**2  +BF  **2 

Fl  =  SQRTF ( YY  >  *COSF ( YF ) 

F2=-SQRTF(YY)*SINF(YF) 

P2«B ( 1 ) *OMEGA 


oo 


Ql  =  -D( 1 )  *OMEGA**2+G  ( 1 ) 

Q2=E ( 1 ) *0MEGA 

Rl=-D ( 2 ) *OMEGA**2+G ( 2 ) 

R2=E ( 2 ) *0MEGA 

Sl  =  -A<2  >*OMEGA**2+C(2) 

S2=B ( 2 ) *0MEGA 

SMQ 1=SM1*Q1-SM2*Q2 

SMQ2=SM1*Q2+SM2*Q1 

FS1=F1*S1-F2*S2 

FS2=F1*S2+F2*S1 

QRl=Ql*Rl-Q2*R2 

QR2=Q1*R2+Q2*R1 

PSl=Pl*Sl-P2*S2 

PS2=P1*S2+P2*S1 

FRl=Fl*Rl-F2*R2 

FR2=F1*R2+F2*R1 

SMP1=SM1*P1-SM2*P2 

SMP2=SM1*P2+SM2*P1 

DEN  =(  QR 1 -PS I T** 2'+ ruR2-FS 2  )  **2  -  . . •-  ' 

ZZl= ( (SMQ1-FS1 )*(QR1-PS1 )  +  ( SMQ2-FS2 )*(QR2-PS2  )  )  /DEN 
ZZ2= ( (SMQ2-FS2)*(QR1-PS1 )-( SMQ1-FS1 ) * ( QR2-PS2 ) l/DEN 
DENG- < i QR1-PS1 > **2+(QR2-PS2  )  **2 ) 

SMM 1  =  (  (  FR1-SMP1 )*<QR1-PS1 )  +  ( FR2-SMP2 ) * ( QR2-PS2 ) )/DENO 
SMM2=( ( FR2-SMP2 ) * ( OR 1-PS1 ! - < FR1-SMP1 ) * ( QR2-PS2  ) )/DENO 
Z=5QRTF(ZZ1**2+ZZ2**2) 

THFT  =  SQRTF<  SMM1**2  +  SMM2**2  ) 

FLZ  =  ATANF ( ZZ2/ZZ1  ) 

FLTH=ATANF(SMM2/SMM1) 

PRINT  7, ( V ( I ) .Z.THET.FLZ.FLTH) 

7  FORMAT  (5F15.8) 

9  1=1+1 

I F ( 1-18)53.53.19 
19  CONTINUE 
CALL  EXIT 
END 
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*HEAVE  AND  PITCH  WITHOUT  CONTROLS  LAMDA  =  4F T * T I  ME  HISTORY  SOLUTION 

*  XEQ 

*  L I S  T  8 

C  HEAVING  AND  PITCHING  RESPONSES 

DIMENSION  V(30) >A1 (30) *A2( 30) *B1 (30) *B2( 30) »A(2>  *B<2) *C(  3 ) >0(2 ) * 
2E ( 2 ) *G ( 2 ) 

READ  1» ( V( I ) »B1 ! I ) »B2( I ) »A1 (  I ) *A2( I ) *1  =  1 »9> 

1  FORMAT  (5F14.8) 

J=1 

1=1 

C ( 1 > =7.74*1.99 
C ( 2 ) =2 . 78 
G( 2 ) =2.43 
D ( 2 )  =0.0251*1.99 
D ( 1 (=0.0251*1.99 

2  V=V(I) 

20  A ( 1 ) =A1 ( I ) 

A ( 2  >  =A2  <  I  ) 

33  B  ( 1  ) =B 1  (  I  ) 

B ( 2 >=B2 (I) 

C ( 3 ) =-0 . 4*V ( I ) +0. 05*V ( I ) **2 
E ( 1 ) =0 . 2*V ( I ) -0 .0 1*V ( I ) **2 
E ( 2 )=-0.1*V( I ) 

G ( 1 >  =0.076*1. 99*V( I ) **2 
OMEGA= 1 . 5  7* ( V ( I  1-4.54) 

70  AF  =  0. 072-0  ,s  0081 6*  V(  I) 

BF=0. 13+0. 04*V ( I ) 

CM=-0.9+0.07*V( I ) 

DM=-0. 042+0. 0008*V( I ) 

8  CALL  RESPON  (Z*THET*A*B»C*D*E*G*OMEGA*AF*BF*  CM  » DM  *  V ) 

13  IF(J-1 ) 14*16.14 

16  E ( 1 ) = — 0 • 3  * V ( I ) +0. 0 1*V ( I )**2 
E  <  2 ) =0 . 1*V ( I ) 

J  =  0 

GO  TO  8 

14  1=1+1 
J  =  1 

I F ( 1-912.2.3 

3  CONTINUE 
CALL  EXIT 

end 
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*  LISTS 

subroutine  RESPONIZ ,THET ♦  A.B*C»D*E.»G»OMEGA»AF »6F.CM*DM»V> 
DIMENSION  A ( 2 ) >  B  (  2  ) * C ( 3 ) » D  S  2  )  *  E  (  2  ) *G(2) 

PRINT  l.V 

1.  FORMAT  !1H  6E20.8) 

T  =  0  •  ' 

Z=0.001 
DZ  =  0 . 

THET=0. 

dthet=o. 

TPRINT=0.0 

7  F=AF*COSF(OMEGA*T >  +BF*SINF(OMEGA*T) 

XM=CM*COSF (OMEGA*T )  +DM*S INF (OMEGA*T ) 

SENSE  LIGHT  1 

FTHET  =  THET/SQRTF<  ABSFSTHET ) ) 

2  DDZ=(A<2)*F-D(1)*XM+(D(1)*E(2)-B(1)*A(2))*DZ  +(D(1>*G(2)  -C ( 1 ) * A 

2(2) ) *  Z  +  ( D ( 1 ) *8 ( 2  )  -E  < 1 ) *A ( 2 ) }*DTHET  +(D(1)*C(2)  -G(1)*A(2>) 
3*THET+A ( 1 )*C(3)*FTHET)/(A( 1 )*A(2)-D(1)*D(2> ) 

DDTHET  = ( -D ( 2 ) *F  +A(1)*XM  +(.B<1)*D(2>  -A  ( 1 )  *E  (  2  )  )  *DZ  +  (  C  ( 1 )  *D  ( 2  ) -A  I 
21)*G(2)*Z+(E(1)*D(2)-A(1)*8(2)*DTHET+(G( 1 ) *D ( 2 ) —A ( 1 )*C<2) )*T 
3HET-A(1)*C(3)*FTHET)/(A(1)*A(2)-D(1)*D(2) ) 

IF  ST  -TPRINT)  4 ♦ 3  f  3 

3  PRINT  1  .T  .Z.THET .DDZ. DDTHET 
TPR INT=TPRINT+0. 1 

4  I F ( T-12. 15,6.6 

5  T=INDVF(T*0.005) 

THET=DPNVF(THET.DTHET) 

Z=DPNVF(Z.DZ) 

DZ=DPNVF( DZtDDZ) 

DTHET=DPNVF(DTHET» DDTHET ) 

GO  TO  7 

6  RETURN 
END 
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*HEAVE  AND  PITCH  WITH  CONTROLS  LAMDA=4FT .LINEARIZED  EQUATIONS 

*  XEQ 

*  L I S  T  8 

*  CONTROL  OF  PITCH  AND  HEAVE 

DIMENSION  V ( 30 ) » A 1 (  30 ) . A2 ( 30 ) »B1 ( 30 ) .82 ( 30 ) . A ( 2 ) »B ( 2 ) >C ( 3 ) *D ( 2  )  . 
2E(  2  >  »G(2  >  .  THETA  (  30)  ,R  (  30  )  .  S  (  30  )  .  T  (  30 )  »U  (  30  ) » DK.1  (  30)  »DK2(  30)  . 

3DK3  <  30  )  » DK4  (  30  )  .PK1( 30)  * PK2  (  30  )  .PK3(30)  . PK4 ( 30 >  >PK1 1 ( 30 )  .PK22130). 
4PK.33  (  30  )  »  PK44  (  30  ) 

PRINT  5 

READ  1. ( V( I ) ,B1 ( I ) »B2< I ) . A1 (I)  .A2  t  I  ) . 1  =  1 .4) 

PRINT  1. ( V< I ) »B1< I ) .62 ( I ) .A1 ( I ) .A2< I ) .1=1.4) 

READ  1. ( THETAt I ). 1=1,4) 

PRINT  1 ♦ ( THET  A ( I ) . 1  =  1 ,4) 

5  FORMAT  (8H  LAMDA=4  ) 

READ  11»(S(I)»T(I).R(I).U( I). 1=1.4) 

READ  111 . (DK1 ( I ) . D<2 ! I ) »D<3( I ) » DK4 ( I ) » PK 1 1 ( I ) .PK22 ( I ) . 1  =  1 .4) 

READ  111 , ( PK33 ( I )  »P<44( I ) .PKl(I) »PK2( I) »PK3< I ) » PK4 ( I ) . 1  =  1 .4) 

1  FORMAT  (5F14.8) 

II  FORMAT  (4F10.5) 

III  FORMAT  (6F10.5) 

PRINT  30 

30  FORMAT  (75H  K1  K2  K3  K4  Z  THE 

2T  A  XI  ) 

1=3 

C( 1 ) =7.74*1 .99 
G ( 2 ) =2 . 3  5 
56  PKK 1 =PK 1 (  I  ) 

PKK  2  =  PK2 ( I  ) 

PKK 3  =  PK3 (  I  ) 

PKK4=PK4( I  ) 

53  PRINT  2  » V (  I  ) 

2  FORMAT  ( 7H  V(I)=F15.8) 

*  WITHOUT  CONTROLS 

C( 3 1=0.09 17*V( I )**2-0.692*V( I ) “0.493 
C  (  2  )  =  (  2  •  78*THET  A  (  I  ) +C  ('3  )  *SQRTF  (  THET  A  (  I  )  )  ) /THETA  (  I  ) 

G ( 1 ) =0.076*1 .99* V ( I ) **2 
OMEGA=l . 57* ( V ( I ) -4 . 54 ) 

70  AF  =  0. 072-0. 0C816*V(  I  ) 

BF  =  0 .13  +  0 . 04*V (  I  ) 

CM=-0.9+0.07*V( I ) 

DM=-0. 042+0. 0008* VC I ) 

20  A ( 1 ) =A1 ( I  ) 

A ( 2  )  =A2 (  I  ) 

B ( 1 ) =B 1 (  I  ) 

B ( 2  )  =B2  t  I  ) 

C(3)=0.0917*V(I ) **2-0 . 692* V (  I  1-0.493 
Dt 1 ) =0.0251*1.99 
DC  2 >=D( 1  ) 

E  1 1  ) =-0.3*V( I )+0.01*V( I )  **2 
E ( 2 >=-0.1*V( I ) 

*  with  foils 

AF= AF+0 • 01 1*V ( I ) 

CM=CM+0 . 022*V ( I ) 

A ( l ) = A ( 1 )+0.012+0.003*V(I ) +0. 103*V ( I ) **2*0 . 1* < PKK2+PKK4) 

A ( 2 ) =A  <  2 ) +0.024+0 .006*V( I ) +0 . 206*V ( I ) **2*0 . 1*PKK1 
B ( 1 )=B( 1 )+0.204*V( I )-0.103*V( I ) **2* ( PKK2+PKK4 ) 
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B<2)=B(2)+0«408#V( I )-0.204*V< I >**2*PKK1 
C ( 2  >  =C ( 2  )  +0.206*V (  I)**2 
D  ( 1 ) =D ( 1 ) +0 « 018 
D ( 2 ) =D ( 1  i 

E( 1 )=E( 1 )+0.204*V ( I )-0.104*V( I ) **2* ( PKK 1  +  PKK.3 ) 

E  <  2 )=E( 2  )+0.204*V< I  1-0.208 *V( I )**2*PKK2 
G(  1  )=.G(  1  )+0.206*V(  1  )  **2 

*  SOLUTION  OF  SET  OF  DIFFERENTIAL  EQUATIONS 

XF=ATANF(DM/CM) 

XX=CM**2+DM**2 

SMl=SQRTF(XX)*COSF(XF) 

SM2=-SQRTF(XX) *SINF(XF) 

YF=ATANF ( BF/AF ) 

YY=AF**2+BF**2 

F1=SQRTF(YY)*C0SF(YF) 

F2=-SQRTF (YY)*SINF(YF) 

*  STABILITY  CRITERIUM 

21  X1=B(1 }*C<2 )+C( 1)*B(2)-E(2)*G(1 )-G(2!*E(l) 

IF ( XI ) 3  *  3  *22 

22  X2  =A  1 1 ) *C  <  2 ) +C  < 1 ) *A  <  2  >  +B  < 1 ! *B ( 2 ) ~D ( 2  >  *G  < 1 ) ~G ( 2 ) *D  < 1 ) -E ( 1 >  *E  <  2 ) 
I F(X2 >3*3*23 

X3=B ( 1 >*A(2)+A( 1>*B(2)-D(2 >*E(1 )-E(2)*D( 1) 

I F(X3 >3*3.24 

24  XO  =  C ( 1 ) *C ( 2  >  — G ( 1 >  *6(2 ) 

X4=A(1)*A(2)-D(1)*D<2) 

X5=X1*X2-X3*X0 
IF(X5 >3*3*25 

25  X6=X1*( X2*X3-X4*X1 )~X3**2*X0 
IF(X6)3.3.26 

26  P 1 = -A ( 1 ) *0MEGA**2+C ( 1 ) 

P2=B(1)*0MEGA 

Ql=-D( 1 ) *0MEGA**2+G ( 1 ) 

Q2=E ( 1 ) *0MEGA 

Rl=-D( 2 )*OMEGA**2+G(2) 

R2=E ( 2 ) *OMEGA 

Sl=-A ( 2 ) *0MEGA**2+C ( 2 ) 

S2  =  B ( 2 ) IOMEGA 

SMQ1=SM1*Q1-SM2*Q2 

SMQ2=SM1*Q2+SM2*Q1 

FSl=Fl*Sl-F2*S2 

FS2=F1*S2+F2*S1 

QRl=Ql*Rl-Q2*R2 

QR2=Q1*R2+Q2*R1 

PS1=P1*S1-P2*S2 

PS2=P1*S2+P2*S1 

FRl=Fl*Rl-F2*R2 

FR2=F1*R2+F2*R1 

SMP1=SM1*P1-SM2*P2  '  •' 

SMP2=SM1*P2+SM2*P1 
DEN=(QR1-PS1)**2+(QR2-PS2! **2 

ZZl=(  (SMQ1-FS1  )*(QRl-PSl>+(SMQ2-FS2)*<QR2-PS2)  )/DEN 
ZZ2=( (SMQ2-FS2)*(QRl-PSl)-<  SMQ1 -FS1 ) * ! QR2-PS2 ) ) /DEN 
D£nO= ( ( QR1-PS1 )**2+(QR2-PS2 )**2 ) 

SMM 1= ( ( FR1-SMP1)*(QR1-PS1>+(FR2~SMP2)*(QR2-PS2 > )/DENO 
SMM2=( ( FR2-SMP2)*<QR1-PS1)-(FR1-SMP1)*(QR2-PS2) > /DEMO 
Z=SQRTF<ZZ1**2+ZZ2**2) 
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THET=SQRTF(SMM1**2+SMM2**2 5 
I F ( Z-S ( I ) 550.50.3 

50  I F { THET-T (I)>51.51.3 

51  PRINT  7  » (  PKK1 » PK.K2  . PKK3 »PKK4 »Z ♦  THET .XI) 

7  FORMAT  (7F10.5) 

I F  C  Z-R (1)58.8.3 
6  IF(THET-U(I) 518.18.3 

3  PK.K  1  =  PKK  1-DK1  (  I  5 

I F ( PKK1-PK11 ( I) 515.15.20 

15  PK<2  =  PKK2-DK2(  I  5 
PKK1  =  PK1 (  I  5 

I F ( PKK2-PK22 ( I  5  516.16.20 

16  PKK3  =  PKK3+DK3 (  I  ) 

PKK 1  =  PK 1 (  I  5 
PKK2=PK2( I  5 

I F ( PKK3-PK33 ( I  5  5  20.1  7.17 

17  PKK4  =  PKK4-DK4(  I  5 
PKK1  =  PK1 (  I  5 

PKK  2  =  PK2 ( I  5 
PK<3  =  PK3(  I  5 

I F ( PKK4-PK44 ( I  5  5 18.18,20 

18  1=1+1 

119  IF ( 1-4) 56 ,56* 19 

19  CONTINUE 

call  EXIT 

end 


*H£AVE  And  PITCH  WITH  CONTROLS  LAMDA=8FT*TIME  history  solution 

*  XEQ 

*  LIST8 

C  HEAVING  AND  PITCHING  RESPONSES 

DIMENSION  V(30 ) *A1 ( 30) *A2( 30) »B1 ( 30) *B2 (30) » A ( 2 ) * B ( 2 ) »C(3 ) »D(2 ) » 
2E( 2 ) >G< 2 ) *PK1 ( 30) *PK2< 30) .  PK31 30) *PK4< 30) 

READ  1, <V( I ) .61 ( I ) .82 ( I ).A1 t I ) »A2( I ) .1=1.12) 

READ  11* ( P  K 1 ( I ) » P  K  2 ( I ) *  PK3  ( I ) *  PK4 ( I ) » 1  =  1 *12 ) 

11  FORMAT  (4F10.5) 

1  FORMAT  (5F14.8) 

PRINT  20. (  I  *A1 ( I )  .  1=7.12) 

20  FORMAT  (4H  A1 (  I  2 . 2H ) =F10. 5  1 
C ( 1 ) =7.74*1.99 
1=7 

*  WITHOUT  CONTROLS 
C(  2  1=2.78. 

D ( 1 )  =0.0251*1.99 
D !  2- 1  =  D  ( 1 1 
G<  2 ) =2.43 
PRINT  15 

15  FORMAT  ( 8H  LAMDA=8 1 

2  V=V  <  I  ) 

56  PKK1=PK1(I1 

PK<2=PK2 (  I  1 
PK<3=PK3<  I  ) 

PK<4  =  PK4 ( I  ) 

A ( 1  )  =A1 (  I  ) 

A ( 2  )=A2 (  I  ) 

B ( 1 )  =  ( 0 . 01 5  5*V ( I ) **2-0 . 083*V ( I  1+0.152  1*1.99 

B( 2 )  =  (-0.000147*V( I )**3+0.01603*V( I  1 **2-0. 0372*V ( 11+0.02471*1.99 
7  El  1 )=-0.5*V( I )+0.02*V( I  1**2 

EC  2)=0.1*V( 1 1 
Gill =0. 076*1. 99*V< 11**2 
OMFGA=1.57*(V( 11-6.41/2. 

AF=  0 • 12-0 . 01*V ( I  1 
BF  =  -1 . 1 +0 . 1 5*V (  I  1 
CM  =  ~1 • 5-0 . 1 5*V ( 1 1 
DM=0. 094-0. 0206*V( I  1 

*  WITH  FOILS 

AF  =  AF-0.041*V( I  1 

DM=DM+0 . 124*V ( I  1 

All )=A( 1 1 +0.01 9+0. 004*V( I  1 

A(2l=A(2)+0. 038+0. 008 *Vt I) 

B ( 1 1 =B ( 1  1 +0.292*V ( I  1 
B(2)=B(2)+0.68*V( I ) 

C ( 3 1 =  — 0 . 4*V ( I  1+0.0  5*V ( I  1**2 

G( 1 )=G( 1  )+0.28*V( I  1**2 

Cl  2  )=C(  2  )-0.336*V(  I  1  ** 2*P<K.2 

13  CALL  RESPO  (2*THET.  ,A».B.*C*D*E,G*  OMEGA  *AF»  BF  *  CM ,DM  > V  *PKK1 ,PKK2  *PKK3 
2.PKK41 
6  1=1+1 

I F ( 1-12)2*2,3 

3  CONTINUE 
CALL  EXIT 
END. 
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*  LISTS 

SUBROUTINE  RESPO  (Z.THET.A.B.C.D.E.G.OMEGA.AF.BF.CM.DM.V.PKKl ,pkk 
22.PKK3.PKK4) 

DIMENSION  A ( 2 )  * B ( 2 ) *C(3)»D(2)*E(2) * G ( 2 ) 


PRINT 

1  .V 

1 

FORMAT 

(1H  F20.8) 

PRINT 

12 

12 

FORMAT 

( 1 1  OH  T 

Z 

THET  DELP 

2LH 

DZ  DTHET 

DDZ 

DDTHET 

T  =  0. 

Z=0.001 
D2  =  0. 

THET=0. 

dthet=o. 

tprint^o.o 

7  F=AF*COSF(OMEGA*T)  +BF*SINF(OMEGA*T ) 

XM=CM*COSF< OMEGA* T )  +DM*SI NF ( OMEGA*T ) 

SENSE  LIGHT  1 

F=F+0.12*V**2*PKK4*DZ/ABSF(DZ) 

XM=XM+0.336*V**2*PKK1*DTHET/ABSF (DTHET ) 

FTHET  =  THET/SQRTF( ABSF(THET)  ) 

2  DDZ= ( A ( 2 ) *F-D< 1 )*XM+( D( 1 >*E< 2 ) -B  C 1 >*A(2 )  )*DZ  +(D(1)*G(2>  — C ( 1 )*A 
2(2)  )  *Z  +  (  D  ( 1 )  *B  (  2  )  -E ( 1 )  #A ( 2  )  ) *DTHET  +(D(1)*C(2>  -G ( 1 ) *A ( 2  )  ) 
3*ThET+A( 1 )*C(3)*FTHET) / ( A ( 1 ) *A ! 2 ) -D C 1 ) *D ( 2 ) ) 

DDTHET  = ( -D ( 2 ) *F  +A(1)*XM  +  (  B  <  1 )  *  D  (  2  )  -A ( 1 ) *E ( 2 ) ) *DZ+ ( C ( 1 ) *D ( 2 ) -A ( 
21)*G(2)*Z+(E( 1 )#D(2)-A(l)*B(2)*DTHET+(6( 1 )  *D  ( 2 ) “A ( 1 )  *C  ( 2 ) >*T 
3HET-A(1 ) *C ( 3 ) * FTHET ) / <  A ( 1 ) ( 2 ) -D ( 1 ) *D ( 2 ) ) 
DElH=PKK4*DZ/A8SF(DZ)-(0.37/V)*CQSF(OMEGA*T)-THET 
DELP=PKK1*DTHET/ABSF( DTHET )+ (0.37/ V)*S INF ( OMEGA*T ) +THET* ( PKK2-1. ) 
IF  (T  -TPRINT)  4.3*3 

3  PRINT  11 .T.Z.THET. DELP.DELH.DZ. DTHET .DDZ, DDTHET 

11  FORMAT  ( 1H  9F12.5  ) 

TPR7NT=TPRINT+0.1 

4  IF ( T— 10« ) 5.6.6 

5  T= I NDVF ( T .0 • 00  5 ) 

THET=DPNVF(THET, DTHET) 

Z=DPNVF(Z.DZ) 

DZ=DPNVF ( DZ.DDZ ) 

DTHET=DPNVF( DTHET. DDTHET) 

GO  TO  7 

6  RETURN 
END 


u_ 
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FiG.  S'  (ft) 

PITCB2BG  AMPLITUDES  WITH  CONTROLS 
X  =  8  ft 

(The  numbers  1,  2,  3,  b  refer  to  the  same  control  systems  as  In  Pig. 
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